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ABSTRACT 

Aims. The paper is devoted to the study of the underlying stellar population of a sample of 65 nearby early-type galaxies, 
predominantly located in low density environments, a large fraction of which show emission lines. 

Methods. Ages, metallicities, and [a/Fe] ratios have been derived through the comparison of Lick indices measured at different 
galacto-centric distances (7 apertures and 4 gradients) with new simple stellar population (SSP) models that account for the 
presence of a/Fe-enhancement. The SSPs cover a wide range of ages (10^ - 16 X 10** yr), metallicities (0.0004 < Z <0.05), and 
[a/Fe] ratios (0-0.8). To derive the stellar population parameters, we use an algorithm that provides, together with the most 
likely solution in the (age, Z, [a/Fe]) space, the probability density function along the age-metallicity degeneracy. 
Results. We derive a large spread in age, with SSP-equivalent ages ranging from a few to 15 Gyrs. Age does not show any 
significant trend with central velocity dispersion ctc, but E galaxies appear on average older than lenticulars. On the contrary, 
a significant increasing trend of metallicity and [a/Fe] with ctc is observed, testifying that the chemical enrichment was more 
efficient and the duration of the star formation shorter in more massive galaxies. These latter two relations do not depend on 
galaxy morphological type. 

We have also sought possible correlations with the local galaxy density Pxyz, but neither metallicity nor a-enhancement show 
clear trends. However, we find that while low density environments (LDE) {pxyz < 0.4) contain very young objects (from 1 Gyr 
to 4 Gyr), none of the galaxies in the higher density environments (HDE) (40 % of galaxies with a measured density) is younger 
than 5 Gyrs. Considering the lack of environmental effect on the [a/Fe] ratio and the high value of [a/Fe] in some young massive 
objects, we argue that young galaxies in LDE are more likely due to recent rejuvenation episodes. By comparing the number 
of "rejuvenated" objects with the total number of galaxies in our sample, and by means of simple two-SSP component models, 
we estimate that, on average, the rejuvenation episodes do not involve more than 25 % of the total galaxy mass. 
The good quality of the data also allow us to analyze the gradients of stellar populations. An average negative metallicity 
gradient AlogZ/Alog(r/re) ~ —0.21 is firmly detected, while the age and a-enhancement spatial distributions within re/2 
appear quite flat. These observations suggest that, in a given galaxy, the star formation proceeded on similar timescales all 
across the central re/2 region, but with an increasing efficiency toward the center. 



Key words. Galaxies: elliptical and lenticular, cD - Galaxies: fundamental parameters - Galaxies: formation - Galaxies: evolution 



1. Introduction 

This is the third paper of a series dedicated to the study 
of early-type galaxies with emission lines (Rampazzo et 
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al. 120051 hereafter Paper I; Annibali et al. 120051 here- 
after Paper II). A recent work by Falcon-Barroso et al. 
([20061) (see also Sarzi et al. [2006|) based on SAURON ob- 
servations shows that the incidence of ionized-gas emis- 
sion in early- type galaxies is about 75%. By morphological 
type, lenticular galaxies display a slightly higher content 
of ionized gas than elliptical galaxies (83% versus 66%), 
in good agreement with the pioneering survey performed 
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by Macchetto et al. I|1996p . who found emission in 85% 
of the lenticulars and 68% of the elhpticals in their sam- 
ple. Falcon-Barroso et al. (J20065 found similar percentages 
when their sample was divided according to environment 
(83% field, 66% cluster). Furthermore, the incidence of 
ionized gas shows no correlation with either luminosity or 
the presence of a bar in the galaxy. The above percentages 
suggest that the phenomenon of emission lines is quite 
common among early-type galaxies: early-type galaxies 
without emission lines are more exceptions rather than 
the rule in this class of galaxies. The first step to under- 
standing the evolution of these galaxies and the nature of 
their ionized gas consists in studying the galaxy underly- 
ing stellar population. 

For this purpose we make use of the of absorption- 
line indices that have been proven to be a power- 
ful tool to disentangle age and metallicity effects (e.g., 
Buzzoni et al. [T992l Worthey [1992'; Gonzalez [1993] (G93); 
Worthey \WM Buzzoni et al. 11994; Worthey et al. [TOMl 
Bressan et al. [TOMl Leonardi & Rose \WM Worthey & 
Ottaviani lTMTl Trager et al.lME; Longhetti et al. HMS} 
Vazdekis[inMl Longhetti et al. J999; Longhetti et allWUl 
Trager et al. I2000al Kuntschner [20001 Beuing et al. 2002; 
Kuntschner et al. [20021 Mehlert et al. [2003)) . More recently 
the interpretation has been faced with the more difhcult 
task of also accounting for the presence of non-solar scaled 
abundance patterns in early- type galaxies, and of deriving 
age, metallicity, and [a/Fe] ratio at the same time (e.g., 
Tantalo et al. [T9961 Tantalo & Chiosi I2004al Thomas et 
al. 120031 Thomas et al. I2005p . Indeed, it has been shown 
that the a/Fe element ratio plays a fundamental role in 
deriving the formation epochs. This ratio, in fact, quan- 
tifying the relative importance of Type II and Type la 
supernovae in the enrichment of the ISM, carries infor- 
mation about the timescale over which the star formation 
occurred (Greggio & Renzini 1983', Matteucci & Greggio 
UW6{ Pagel & Tautvaisiene IT995. Thomas et al. [TgM|) . 

In this paper we try to constrain the epoch of the 
baryon assembly of our sample of 65 nearby early-type 
galaxies predominantly located in low density environ- 
ments, a large fraction of which show emission lines (Paper 
I and II), and seek possible correlations with their mor- 
phological type (elliptical and lenticular) , velocity disper- 
sion, and environmental density. A large part of this paper 
is dedicated to the development of our own a-enhanced 
models and to the comparison with the models already 
presented in the literature. The stellar population param- 
eters are derived by means of a new algorithm based on 
the probability density function. The algorithm provides, 
together with the most probable solution in the (age, Z, 
[a/Fe]) space, the solutions along the age- metallicity de- 
generacy that are within 1 a error from the observed index 
values. Our analysis is based on the use of three different 
line-strength index diagnostics, including the Mgb index, 
the <Fe> index, and a Balmer line (II/3, II7, or US). The 
use of higher order Balmer lines allows us to minimize 
spurious effects that may arise from uncertainties in the 
emission correction of the H/3 absorption line. The pa- 



per is organized as follows. Section 2 presents an overview 
of the sample, summarizing its main properties and fo- 
cusing on the delicate step of the correction for emission 
infilling of the most commonly used age indicator, the II/3 
absorption line. In Sect. 3 a complete description of the 
procedure adopted to compute the new a-enhanced SSPs 
is provided. Section 4 describes how the new SSP models 
have been implemented within an algorithm that allows 
us to derive the stellar population parameters, i.e., age, 
metallicity, and a/Fe ratio. Section 5 presents the results. 
In particular we analyze both the stellar population pa- 
rameters derived for the apertures and gradients, and their 
relations with velocity dispersion and local galaxy density. 
In Sect. 6 we discuss our results and compare them with 
the literature. Finally, a summary and our conclusions are 
presented in Sect. 7. 

2. Further analysis of the line-strength indices 
data set 

2.1. A sample overview 

The procedure of data reduction and analysis of the 65 
nearby early-type galaxies studied here has been already 
discussed in Papers I and II. Table 1 summarizes the main 
characteristics of the sample. Column (1) gives the galaxy 
identification name; Col. (2) and (3) provide the galaxy 
morphological classification according to RSA (Sandage 
& Tammann [TM7[) and RC3 (de Vaucouleurs et al. [TMT)) . 
respectively: only in few cases do the two catalogues dis- 
agree in the distinction between E and SO classes; Col. (4) 
gives the galaxy systemic velocity, Yhei, which is lower 
than ~5000 km s~^; Col. (5) provides the richness param- 
eter pxyz fTullv [TO88p : it represents the density of galaxies 
brighter than -16 apparent B-mag in the vicinity of the 
entry, in galaxies xMpc~^. The local density was deter- 
mined on a three-dimensional grid with 0.5 Mpc spacings. 
A Gaussian smoothing function was used to calculate the 
contribution of each member of the sample brighter than 
-16 B-mag to the density of the specific location. In the 
TuUy catalogue, the galaxies (considered in the local envi- 
ronment) within 40 Mpc (3000 km s~^) with an apparent 
magnitude B-mag <-16 arc 2189. If one entry is at a dis- 
tance larger than 40 Mpc, the local density is not recorded. 
The galaxies of our sample are mainly located in low den- 
sity environments. The local density of our galaxies varies 
from pxyz ~ 0.1, characteristic of very isolated galaxies, to 
Pxyz ~ 4, which is characteristic of denser galaxy regions 
in the Virgo cluster. For comparison, in the TuUy 119881 
catalogue, objects like NGC 1399 and NGC 1389, Fornax 
cluster members, have values of pxyz = ^.5d and 1.50, re- 
spectively. Thus, our sample, even though biased towards 
low density environments, contains a fraction of galaxies 
in relatively dense environments. 

The sample spans a large range in central velocity dis- 
persion, « from 115 to 340 km s~^ (see Papers I and II 
for detailed descriptions). Our long-slit spectra cover the 
(3700 - 7250) A wavelength range with a spectral resolu- 
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tion of «7.6 A at 5550 A. For the detailed description of 
the data reduction procedure and of the derivation of line- 
strength indices in the Lick system, we refer to Papers I 
and II. We only recall here that for each galaxy, 25 Lick 
indices (21 belonging to the original Lick set (see Trager et 
al. (|1998p for the pass-band definitions) plus 4 higher order 
Balmer lines introduced by Worthey & Ottaviani (|1997p ) 
have been measured for 7 luminosity weighted apertures 
(with radii: L5", 2.5", 10", re/10, re/8, re/4, and re/2), 
corrected for the galaxy ellipticity, and 4 gradients (0 < 
r <re/16, re/16 < r <re/8, re/8 < r <re/4, and re/4 < r 

<l-e/2). 

2.2. Hf3 correction for emission infilling 

Because of the presence of emission in a significant frac- 
tion of the sample, the problem of possible contamination 
of the H/3 absorption line deserves particular attention. In 
Papers I and II we have shown two different methods to 
correct the H/3 index for possible emission infilling, based 
respectively on the measure of the [OIII]A5007 A and of 
the Ha emission lines. The first method is based on the 
correlation found by G93 between the B./3 and the [OIII] 
emission in his sample of early-type galaxies, such that 
EW(H/3om)/EW([OIII]A5007) = 0.7. The second method 
rests on the possibility of measuring the Ha emission at 
6563 A, from which the emitted flux in the H/3 line results 
in Fh/3 = l/2.86FHa (Osterbrock 198¥). In Paper I and II, 
we have shown that the corrections derived from the two 
different methods are statistically similar. We then pro- 
vided the final H/3 indices, applying the correction from 
the [OIIIJA5007. We noticed, however, that the scatter of 
the relation derived by G93 is large, and that small vari- 
ations in the H/3 index correspond to large differences in 
the derived stellar population parameters when compared 
with models. 

A comparison of the emission corrections derived from 
the [OIIIJA5007 and Ha lines, respectively, is also per- 
formed by Denicolo et al. (2005a) for their sample of early- 
type galaxies. The authors show that the differences in the 
emission corrections derived through the two methods are 
small for the majority of galaxies, but argue that it is 
very dangerous to rely on the use of the [OIIIJA5007 alone 
because the [OIIIJA5007/H/3 ratio may be very different 
from galaxy to galaxy. They use the Ha to correct the H/3 
index. 

The uncertainties concerning the emission corrections 
derived from the Ha line mainly rest on the assumption of 
a value for the Ha in absorption. Denicolo et al. (|2005aP 
adopt an absorption value on the basis of the Ha EW mea- 
sured for their sample of stars (mainly K). Differently, as 
described in Papers I and II, we computed the emitted 
flux in the Ha line using, as a reference, the spectrum of 
an elliptical galaxy (NGC1426) lacking emission lines or 
dust. This galaxy is located in the low tail of the Mg2-(T 
relation. However, we verified that the use of a giant el- 
liptical reflects in a difference of 0.03 A in the value of the 
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R,C3 




Pxyz 








km s ^ 


Gal. Mpc-^ 



computed H/3 emission. Denicolo et al. p005ap allowed 
for a reasonable range of Ha absorption strengths and de- 
termined that the difference in the emission correction is 
smaller than 0.1 A. In this paper we adopt the direct esti- 
mate derived from the Ha emission line for the correction 
of the H/3 index. The adopted H/3 values can be easily re- 
covered from Papers I and II, where we provided, together 
with the final indices in the Lick system, the H/3 emission 
estimates derived both from the [OIII]A5007 and Ha lines. 
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S02{8) poc 


SO pec sp 


4227 




NOG 777 


El 


El 


5040 




NGG 1052 


E3/S0 


E4 


1475 


0.49 


NGG 1209 


E6 


E6: 


2619 


0.13 


NGG 1297 


S02/3{0) 


SABO pec: 


1550 


0.71 


NGG 1366 


E7/S01(7) 


SO sp 


1310 


0.16 


NGG 1380 


S03{7)/Sa 


SAO 


1844 


1.54 


NGG 1389 


S01{5)/SBD1 


SAB(s)0-: 
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1.50 


NGG 1407 


E0/S01(0) 


EO 


1766 


0.42 


NGG 1426 


E4 


E4 


1443 


0.66 


NGG 1453 


EO 


E2 


3906 




NGG 1521 


E3 


E3 


4165 




NGG 1533 


SB02{2)/SBa 


SBO- 
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0.89 


NGG 1553 


S01/2(5)pec 


SA(r)0 


1280 


0.97 


NGG 1947 


S03{D) poc 


SO- poc 


1100 


0.24 


NGG 2749 


E3 


E3 


4180 




NGG 2911 


SOp or S03{2) 


SA(s)0: pec 


3131 




NGG 2962 


RSB02/Sa 


RSAB{rs)0-|- 


2117 


0.15 


NGG 2974 


E4 


E4 


1890 


0.26 


NGG 3136 


E4 


E: 


1731 


0.11 


NGG 3258 


El 


El 


2778 


0.72 


NGG 3268 


E2 


E2 


2818 


0.69 


NGG 3489 


S03/Sa 


SAB{rs) + 


693 


0.39 


NGG 3557 


E3 


E3 


3038 


0.28 


NGG 3607 


S03{3) 


SA(s)0: 


934 


0.34 


NGG 3818 


E5 


E5 


1701 


0.20 


NGG 3962 


El 


El 


1822 


0.32 


NGG 4374 


El 


El 


1060 


3.99 


NGG 4552 


S01{0) 


E 


322 


2.97 


NGG 4636 


E0/S01(6) 


EO-1 


937 


1.33 


NGG 4696 


{E3) 


E + 1 pec 


2958 


0.00 


NGG 4697 


E6 


E6 


1241 


0.60 


NGG 5011 


E2 


El-2 


3104 


0.27 


NGG 5044 


EO 


BO 


2704 


0.38 


NGG 5077 


S01/2{4) 


E3+ 


2764 


0.23 


NGG 5090 


E2 


E2 


3421 




NGG 5193 


S01{0) 


E pec 


3711 




NGG 5266 


S03{5) poc 


SAG-: 


3074 


0.35 


NGG 5328 


E4 


El: 


4671 




NGG 5363 


[S03{5)] 


10: 


1138 


0.28 


NGG 5638 


El 


El 


1676 


0.79 


NGG 5812 


EO 


EO 


1930 


0.19 


NGG 5813 


El 


Bl-2 


1972 


0.88 


NGG 5831 


E4 


B3 


1656 


0.83 


NGG 5846 


S01{0) 


E0+ 


1709 


0.84 


NGG 5898 


S02/3{0) 


BO 


2267 


0.23 


NGG 6721 


El 


E+: 


4416 




NGG 6758 


E2 (morger) 


E + : 


3404 




NGG 6776 


El pec 


E + poc 


5480 




NGG 6868 


E3/S02/3(3) 


E2 


2854 


0.47 


NGG 6875 


SO/a{morger) 


SAB(s)0- pec: 


3121 




NGG 6876 


E3 


E3 


3836 




NGG 6958 


R?S01(3) 


B-l- 


2652 


0.12 


NGG 7007 


SD2/3/a 


SAO-: 


2954 


0.14 


NGG 7079 


SBa 


SB(s)0 


2670 


0.19 


NGG 7097 


E4 


E5 


2404 


0.26 


NGG 7135 


SOI pec 


SAO- pec 


2718 


0.32 


NGG 7192 


S02{0) 


E + : 


2904 


0.28 


NGG 7332 


S02/3{8) 


SO pec sp 


1207 


0.12 


NGG 7377 


S02/3/Sa pec 


SA{s)0-|- 


3291 




IG 1459 


E4 


E 


1659 


0.28 


IG 2006 


El 


E 


1350 


0.12 


IG 3370 


E2 pec 


E2-I- 


2934 


0.20 


IG 4296 


EO 


E 


3762 




IG 5063 


SD3{3)poc/Sa 


SA(s)0+: 


3402 
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3. Modelling line-strength indices for simple 
stellar populations 

3.1. Standard models 

Following the procedure described in Bressan et al. (|1996p , 
to which we refer for details, we have derived line-strength 
indices for SSPs. Indices are constructed by means of a 
central band-pass and two pseudo-continuum band-passes 
on either side of the central band. Molecular bands are ex- 
pressed in magnitude, while atomic features are expressed 
in equivalent width (EW). The definition in EW is: 

Iew = {1-Fr/Fc)AX, (1) 

while the definition in magnitude is 

Imag^-2.5l0g{FR/Fc), (2) 

where Fr and Fq are the fluxes in the line and in the 
pseudo-continuum, respectively, and AA is the width of 
the central band. The flux Fc is obtained by interpolating 
the fluxes in the blue and red pseudo-continua bracketing 
the line of interest to the central wavelength of the absorp- 
tion band. Wc used the passband definitions provided in 
Trager et al. (|1998| for the original 21 Lick indices, and in 
Worthey & Ottaviani (|1997p for the higher-order Balmer 
lines. The integrated indices for SSPs are based on the 
Padova library of stellar models (Bressan et al. I1994p and 
accompanying isochrones (Bertelli et al. I1994p . The SSPs 
indices are calculated using the following method. We de- 
rive both the pseudo-continuum flux Fc from the library 
of stellar spectra used by Bressan et al. (|1994p and imple- 
mented by the revision of Tantalo et al. (|1996p . and the 
line strength index Iew or Imag from the fitting functions 
(FFs) of Worthey et al. p994| and Worthey & Ottaviani 
(|1997p . for each elementary bin Alog-^ and AlogTe// 
of a given isochrone in the HR diagram. The residual flux 
in the central feature Fu is recovered from Fc and Iew 
(Imag) by inverting Eqs. (P) and For an SSP of age T 
and metallicity Z, the fluxes in the line and in the contin- 
uum, Fii^ssp{T, Z) and Fc,ssp{T, Z)^ are then obtained 
by integrating Fc and Fr along the isochrone: 

Fk,ssp{T,Z)= cj){M)Fk{M,T,Z)dM, (3) 

J Ml 

where k denotes the flux in the continuum (C) or in the 
fine (R), (/)(M) is the Initial Mass Function (IMF) dcflned 
between and M[/, and Fi.{M,T, Z) is the flux for a 
star of mass M, age T, and metallicity Z. Once the in- 
tegrated fluxes, Ffi ssp and Fc^ssp, are computed, the 
index deflnition (Eqs. (P) and ([2])) is applied to get the 
integrated SSP index back. Our SSP models have been 
computed adopting a Salpeter (|1955p IMF slope between 
0.15 and 120 Mp,. 



3.2. a/Fe-enhanced models 

To compute the SSP models presented in Sect. 3.1, we 
have adopted stellar isochrones based on solar-scaled 
abundances and FFs calibrated on Milky Way stars. For 
these reasons the standard SSPs reflect the solar element 
abundance pattern. However, there are several evidences 
for supersolar [a /Fe] ratios in early-type galaxies (Peletier 
[TOggi Worthey et al. [TWl Davies et al. ITM51 Carollo & 
DanzigerlUMl Bender & Paquet [19951 Fisher et aL [T995l 
Mehlert et al. UMEi J0rgensen [TMSl Kuntschner [MM 
Longhetti et al. 120001 and others) . 

The major effect of non-solar abundance patterns is 
on the stellar atmospheres. For a given effective tem- 
perature, gravity, and total metallicity, the stellar spec- 
trum, and thus the line-strength indices that are mea- 
sured in it, depend on the specific ratios between the 
element abundances. To account for this effect, we ap- 
ply a correction to the index derived from the FFs us- 
ing tabulated index responses to element abundance vari- 
ations. The responses are derived from model atmospheres 
and synthetic spectra (Tripicco & Bell 119951 (hereafter 
TB95); Korn et al. [MJ5I (hereafter K05), Munari et al. 
[20051 (hereafter Mu05)). A second, less important effect is 
on the evolution of the star and on the stellar opacities. 
A fully self-consistent a-enhanced SSP model should in 
principle use a/Fe-enhanced stellar evolutionary tracks 
(Weiss et al. [TM51 Salaris & Weiss HMEi Van den Berg 
et al. Salasnich et al. In the literature, sev- 

eral a-enhanced SSP models have been presented. In all 
models a correction is applied to the index value to ac- 
count for the impact of the a/Fe-enhancement on the 
stellar atmosphere (Weiss et al. I1995| Tantalo, Chiosi & 
Bressan [TM51 Trager et al. I2000al Thomas et al. [MI5I 
(hereafter TMB03), Thomas et ah [20041 Tantalo & Chiosi 
I2004bl Tantalo et al. [20041 Korn et al. [20051) . SSPs adopt- 
ing a/Fe-enhanced stellar evolutionary tracks have been 
produced as well (Thomas & Maraston ,2003; Tantalo & 
Chiosi [^04ap . but current enhanced models seem to over- 
estimate the blueing of the stellar evolutionary tracks. In 
this section we present new SSP models that are based 
on solar-scaled isochrones (Bertelli et al. I1994p and where 
only the direct effect of the enhancement on the stellar 
atmosphere is taken into account. The models are avail- 
able to the pubhc at http://www.inaoep.mx/~abressan 
or www.stsci.edu/~annibali/. 

3.3. a-enhanced chemical compositions 

Following the procedure adopted by TMB03, we assign el- 
ements to three groups: the enhanced group (a-elements), 
the depressed group (Fe-peak elements), and the fixed 
group (elements left unchanged). Taking as a reference 
model a mixture of total metallicity Z and solar parti- 
tions, a new mixture with the same global metallicity Z, 
but with supersolar [a/Fe] ratio is produced by increas- 
ing the abundance of the a-group and decreasing that of 
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the Fe-group in such a way that the two mass fraction 
variations balance. 

If Xa^Q and Xpe^Q are the total mass fractions of the 
two groups for the solar mixture, and X^ and Xp,. are the 
mass fractions in the new mixture, the degree of enhance- 
ment is expressed by the quantity 



[a/Fe] — log log 



Xpe 



Xpe,<i 



(4) 



The quantities Xa.Q and Xpe^Q are computed by summing 
up the mass fractions of the elements within the a and Fe 
groups in the solar-scaled mixture, and thus depend on 
how elements are allocated within the groups. 

For an arbitrary mixture, the i-th element mass frac- 
tion is: 



The solutions fo, and fpo to ^ are: 



fpo 
fa 



Xq, b + Xpcb 



Xa,b 10(["/F<=ln-["/F<=lb) -f Xpcb' 
fp^lO(["/F<=ln-["/Fc]b) 



(8) 



In the case in which the base model corresponds to the 
solar mixture (i.e., X^^b = Xq,^q and Xpo.b = Xpc,©), the 
quantities fa and fpe are given by: 



Xq,0 -I- Xfo,0 



fa = fpolO["/^°l" 



(9) 



Xi — fi 



-X 



i,0i 



(5) 



and the mass fractions in the new mixture are: 



where fi — 1 for solar-scaled abundances, /; > 1 if the 
element is enhanced with respect to the solar-scaled com- 
position, and /i < 1 if the element is depressed. Once the 
element allocation within the three groups and the global 
enhancement [a/Fe] has been defined, we derive the quan- 
tities fi. Since in our assumption all the elements within 
a group are enhanced/depressed by the same amount, it 
is sufficient to compute fa and fpe- 

Let us call Zi, and [a/Fe]b the total metallicity and 
the enhancement of the starting base model, and Z„ and 
[a/Fe\n the same quantities in the new mixture. The en- 
hanced, depressed, and fixed group mass fractions in the 
new mixture (n) are related to the base model (6) quan- 
tities according to: 



Y — f ^" V 

Xpo,n — fpcT^Xpc^bj 
Xo.n = 7^Xo,b 



(6) 



The quantities fa and fpe are derived by solving the set 
of equations: 

Zn = Xc^n -|- Xpo,n + Xo,n = 

= faTJ— Xa,b + fpoTJ— Xpcb + 
■^b 



[a/Fe]„ - [a/Fe], 



"T^^O.b, 
^b 

/ Xq n 



log 
-log 



\Xpo,n 
Xa^b 
Xpc.b 



log 



, / Xa n 



X 



1 

log — 

tPc 



Po,n 



log 
-I- log 
log 



fx. 



VXpc,0 

Xa,0 



Xpc,0 
Xq b 



X 



Pc,b 



(7) 



Xct^n 
Xpo,n 
Xo,ii 



faTJ^Xct,©, 
Z0 

fpe77^Xpo,0, 
Zi0 



z© 



Xo,0 



(10) 



Adopting the solar abundances of Grevesse & Sauval 
(|1998p . and assigning the elements N, O, Ne, Na, Mg, 
Si, S, Ca, and Ti to the enhanced group and Cr, Mn, 
Fe, Co, Ni, Cu, and Zn to the depressed group, we get 
Xa,Q = 1.251 X 10^2 and Xpe,Q = 1.4 x 10"^. 



3.4. Specific response functions 

Specific response functions allow to isolate the effect on 
line-strength indices due to variations of one element at 
once. The first to quantify this effect were TB95 who 
computed model atmospheres and synthetic stellar spec- 
tra along a 5-Gyr-old isochrone with solar metallicity. 
The model atmospheres were computed for three points 
along the isochrone representative of the evolutionary 
phases of dwarfs {Teff — 4575 K, log g — 4.6), turn-off 
(6200 K, 4.1), and giants (4255 K, 1.9). On each syn- 
thetic spectra, TB95 measured the absolute values Iq for 
the original 21 Lick indices. Then they computed new 
model atmospheres, doubling the abundances of the dom- 
inant elements C, N, O, Mg, Fe, Ca, Na, Si, Cr, and Ti 
(A[Xi] = 0.3) in turn, and derived the index changes AJ. 
The index response for an increase of -1-0.3 dex of the i-th 
element abundance is defined as: 



1 dl 

hd\x. 



0.3 



A/' 



TB95 



(11) 



where /q and A/tb95 are, respectively, the index com- 
puted by TB95 for solar composition and the index change 
measured for a -1-0.3 dex abundance increase. Recently 
K05 have extended the work of TB95 to a wide range of 
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Fe4531 


0.600 


0.9362 


TiOl 






C24668 


1.080 


4.1530 


Ti02 


2.722 


0.0092 


H/3 


0.936 


-0.1275 


m 


0.663 


-0.1094 


Fe5015 


0.783 


1.4040 


H7 


0.669 


0.2911 


Mgl 


0.817 


0.0040 


m 


0.461 


0.8014 


Mg2 


0.822 


0.0120 


H7 


0.864 


0.0899 


Mgb 


0.776 


0.6049 









Notes: a and b are the coefficients of the linear transformation 
EWLick = a X EWatmo -l- b adopted to transform the indices 
derived from the model atmospheres of Korn et al. ()2005)) into 
the Lick system. 

metallicities and have included in their analysis the higher 
order Balmer hnes (H7 and US) as well. We have adopted 
these new responses to compute our SSP models with a- 
enhanced composition. 

As the index responses will be used in the following to 
correct solar-scaled SSP models based on FFs, we have 
first of all compared the solar-scaled indices computed by 
K05 for their model atmospheres with the FFs of Worthey 
et al. (|1994p . Wc have limited the comparison to metal- 
licities [Z/H] > 0, as at lower metallicities the FFs, cali- 
brated on Milky Way stars, reflect super solar [a/Fe] ra- 
tios. We observe that while in some cases there is a good 
match between FFs and model atmospheres, in other cases 
large deviations are present (see Fig. 1). 

Mismatches between theoretical models and observa- 
tions may be due both to uncertainties in the input physics 
of the model and to residual effects in the calibration of 
the model into the observed system (Chavez, Malagnini 
& Morossill996|. We then wonder whether model indices 
should be transformed into the Lick system, in analogy 
with the transformation usually performed in the liter- 
ature on measured indices through observation of stan- 
dard stars in common with the Lick library (Worthey 
& Ottaviani I1997p . We derive the linear transformation 
(Iff = a X latmo + b), which is needed to match the at- 
mospheric indices with the FFs. The coefficients a and b 
of the transformation are given in Table 2. In Sect. 3.5 
the index correction will be computed adopting the K05 
responses revised in the following way: 



where Iq is the index from the solar scaled atmospheres, 
la is the index from the a-enhanced atmospheres, and 
I(Q=o,Lick) and AI(c,^Lick) are the solar-scaled index and 
the index change transformed into the Lick system, re- 
spectively. 

An important point to be assessed concerns the metal- 
licity dependence of the indices, and more specifically the 
dependence of the fractional response A///. We recall 
here that the behavior of the equivalent width of a line 
with respect to the abundance of an element is described 
by the curve of growth, and characterized by three dif- 
ferent regimes: a linear part, where EW cx N, for a small 
number of absorbers N (proportional to the metallicity 
Z); a saturated part, where EW cx VlnN, in the Doppler 
regime of saturated lines; and a coUisionally broadened 
regime, described by EW cx -s/N, for large N. 

To assess in which regime the Lick indices are, in a 
first attempt we have analyzed the behavior of the K05 
indices with global metallicity. In principle, such analy- 
sis would require keeping the temperature and the gravity 
fixed and only let the metallicity varying, while the K05 
indices were computed for different (Teg, log g) pairs at 
different metallicities. Nevertheless, T^s and logg varia- 
tions are small within each phase (see Table 2 in K05). By 
analyzing the index behavior, we observe that the Lick 
indices are not linearly related to Z. A viable explanation 
is that, at the Lick-IDS system resolution, the contribu- 
tion to the index absorption strength from weak lines is 
negligible, as weak lines are likely to be fiattened to the 
continuum or blended with stronger lines. In reverse, the 
major contribution to the index equivalent width derives 
from strong saturated absorption lines, for which it would 
be more appropriate to consider the fiat part or the square 
root part of the curve of growth. 

This evidence is strongly supported by preliminary re- 
sults that we have obtained with the last version of the 
ATLAS code (ATLAS12, Kurucz |l993b| , which is based 
on the opacity-sampling (OS) method. Following TB95 
and K05, we have computed model atmospheres and syn- 
thetic spectra for a model with solar abundance ratios and 
for models in which the abundances of the elements C, N, 
O, Mg, and Fe have been changed by different amounts 
([X/H]= -0.2,-h0.2,-h0.3,-h0.5). In Figure [2] we plot the 
index as a function of the logarithmic (left panels) and 
linear (right panels) abundance of the index-dominating 
element for the indices (Mg2, Mgb, Fc5335, Fe5270). For 
more details we refer to Annibali (2005) . We observe that 
the behavior of the indices as a function of the element 
abundances is very well described by a functional form 
that is linear in [X/H] (I = a x [X/H] -I- b) or, equivalently, 
logarithmic in X (I = a x logX/XQ -f- b). 



Al 



(a,Lick) _ (a X la + b) - (a X lo -f b) 



I 



(a=0,Lick) 



a X Iq + b 
a X (Ig - Iq) 
a X In + b 



(12) 



3.5. The Al correction 

As explained in the previous section, specific response 
functions quantify the index change due to the deviation 
of one element from the solar composition. The total index 
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Fig. 1. Comparison between solar-scaled indices obtained by Korn et al. (j2005p for their model atmosphere with 
fitting functions. The comparison is limited to metallicities [Z/H] > 0. The solid line is the one-to-one relation, while 
the dashed line is the performed fit to the data. 



correction for a mixture in which several elements depart 
from the solar abundance path is derived through a linear 
combination of the specific responses, i.e., assuming that 
the total effect on the index can be expressed as the sum 
of the single effects due to the variation of one element at 
once. In the literature, different expressions to derive the 
index change on the basis of the same specific response 
functions have been proposed. 

Using the TB95 responses, Trager et al. (|2000ap pro- 
pose the formula 

f = {n[i+Ro.3(xo]™°-^}-i (13) 

for the fractional index variation, where 
R-o.3(Xi) = gp^.j 0.3 is the fractional index change 
computed by TB95 for an increase of the i-th element 
abundance of 0.3 dex. The computations of TB95 are 
limited to the solar metallicity. Trager et al. (|2000ap 
propose to extend the index correction to non-solar cases 



with the reasonable assumption that the fractional index 
change is independent of metallicity. 

Tantalo & Chiosi (|2004ap adopt the same functional 
form to compute the index variations. TMB03 propose a 
different approach starting from the assumption that the 
line strength index is a linear function of the element abun- 
dance. TMB03 assume that Lick spectral features are in 
the linear regime of the curve of growth. They thus as- 
sume that I oc Xj and consequently Inl cx [Xj] (note, how- 
ever, that the latter implies I oc X"), and then express the 
index through the Taylor series 

Inl = lnI„ + f^|^A[Xi] (14) 
i=i ^^^'J 

However, as discussed in Sect. 3.4, at the Lick resolution 
the major contribution to the index line-strengths arises 
from strong saturated absorption lines, and a logarithmic 
functional form for the index behavior with clement abun- 
dance may be more appropriate than a linear one. 
Then we express the index I through the Taylor expansion: 
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Fig. 2. Mg2, Mgb, Fe5335, and Fe5270 indices as functions of the element abundance. In the left panels the index is 
plotted as a function of the logarithmic element abundance [X/H], while in the right panels as a function of the linear 
abundance X/Xq. Full circles are the model atmospheres computed with ATLAS12, and the solid line connects the 
models, while the dotted line is the performed fit, which is also given in the box. 



I = In 



E 



dl 



5[Xi] 



A[Xi 



The fractional index change is 



AI _ 



dl 

9[Xi] 



A[Xi]. 



(15) 



(16) 



The computations of K05 allow us to derive the frac- 
tional index change ^ at different metallicities. We as- 
sume that we start from a mixture where the total metal- 
licity is Z and the element abundances are 



metallicity Z. If now we perturb the abundance of the i-th 
element, the new mass fraction is given by 



Xi — fiXi_b — fifbTJ — Xi^0, 

and the logarithmic abundance change is 



A[Xi]-log^ = logfi. 



(18) 



(19) 



For a given Z, the expression for the fractional index 
change becomes 



Xi,o — fbTJ— Xi,0, 



(17) 



where fb accounts for the [a/Fe] bias of the stellar libraries 
(see Sect. 3.7), and it is the enhancement/depression fac- 
tor with respect to the solar-scaled abundance for that 



AI, ^Ko.3,z, , 



^ 0.3 

i=l 



(20) 



where Kq.s.z are the new response functions AI/Iq com- 
puted by K05 for an increase of the element abundance of 
-F0.3 dex. 
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Fig. 3. Mgb index versus the other 24 Lick indices for a fixed age of 12 Gyr and metaUicities in the range 0.0004<Z< 
0.05. The metalhcity increases for increasing Mgb values. The long dashed line is the base model. The solid line 
represents the a-enhanced models computed for different [a/Fe] ratios (0, 0.3, 0.5), assuming that the index behavior 
is logarithmic with element abundance. Models with the larger Mgb index correspond to larger [a/Fe] ratios. Full 
squares and open stars are the globular clusters data and the integrated bulge light, respectively, in Puzia et al. 
(|2002p . Open circles are our data sample (Papers I + II). 



3.6. The final SSP model 

With the functional dependence of AI/I on element abun- 
dance described by Eq. ([20|) . the new SSP indices are com- 
puted following the same procedure adopted by TMB03. 

We split the basic SSP model into three evolutionary 
phases, dwarfs (D), turnoff (TO), and giants (G). Gravity 
determines the separation between dwarfs (logg > 3.5) 
and giants (logg < 3.5). Following TMB03, we used the 
fixed temperature of 5000 K to separate the turn-off re- 



gion from the cool dwarfs on the main sequence, indepen- 
dently of age and metallicity. TMB03 have also demon- 
strated that the impact of this choice on the final SSP 
model is small. The Lick indices of the base model are 
computed for each evolutionary phase separately; this is 
done by integrating the continua and line fluxes along each 
phase, and finally deriving the indices Id , Ito , and Iq from 
the continua Fc,d,Fc,tO: and Fc,g and the line fluxes 
Fr,d,Fr^to, and Fr^q of the three phases. At this point 
the index is corrected for each phase, adopting the K05 




Fig. 4. Same as Fig. 3. 

responses in the appropriate phase and computing the AI 
as described in Sect. 3.5. Once the new corrected indices 
Id'ItO' ^^'^ computed, the final a-enhanced SSP 

index is obtained through two different expressions, de- 
pending on if the index is defined as EW or mag. 

EW: 



1 - ^) Fa 



Fc. 



D 



Fc,TO 

X Fc,D + r 



Fee 



AA = 



TO 



C,TO 



+ 1? 



C,G 



Fc,D + Fc,TO + Fc^G 



Magnitude: 



Issp = 



^ _ Fr.,d + Fr^to + Fr^g _ 

Fc,D + Fc,TO + Fc,G_ 



1 - 



Fc\D + Fc,TO + Fc\G 



AA + 



T oci /Fr D + Fr, TO + Fr, g\ 

Issp = -2.5 log — — = 

V * C,D + r c,TO + ^ C,G / 
1 o . 1 lO^oFcD + lO^ToFcTO + lO^^Fc.G 

= 1 —2.5 log ^ — ■ — — — . 

Fc,D + Fc,TO + Fc,G 
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Fig. 5. Same as Fig. 3. 



3.7. The a/Fe bias of stellar libraries 

As discussed by TMB03, an additional complication in 
the construction of the a-enhanced SSP models is due to 
the fact that the fitting functions, on which the standard 
models are based, are constructed from empirical stellar 
libraries that reflect the chemical enrichment history of 
the Milky Way. More specifically, metal poor halo stars 
([Fe/H]< —1), which formed at early epochs in a short 
timescale, are characterized by supersolar [a/Fe] ratios 

0.3). Metal rich disk stars, which formed on longer 
timescales and were enriched by Type la supernovae, in- 
stead exhibit an [a/Fe] trend that decreases from 0.3 to 
solar for increasing metallicity (Edvardsson et al. 119931 
Fuhrmann ll998p . The result is that the base models con- 
structed on the FFs do not have solar abundance ratios 
at every metallicity, but reflect the a/Fe bias of the FFs 
at subsolar metallicity. In Table 3, we show the [a/Fe] 
values adopted for the bias, which are obtained by inter- 
polating the values given in Table 3 of TMB03 to our 
metallicity grid. For a model of given metallicity and to- 
tal [a/Fe], the [a/Fe]b of the bias is the starting point 
for the computation of the fi factors according to Eqs. [8] 
and [5) Summarizing, we have computed new a-enhanced 
SSP models adopting the K05 specific response functions. 
For comparison with TMB03, we present the results for 
models computed for a mixture in which the enhanced el- 
ements are N, O, Ne, Na, Mg, Si, S, Ca, and Ti, while the 
depressed ones are Cr, Mn, Fe, Co, Ni, Cu, and Zn. 

We plot the Mgb index versus the other 24 Lick in- 
dices for a fixed age of 12 Gyr and total metallicities in 
the range 0.0004<Z< 0.05 in Figs. 3, 4, and 5. The metal- 
licity increases for increasing Mgb values. The a-cnhanced 
models represented are computed for different [a/Fe] ra- 
tios (0, 0.3, 0.5). Larger [a/Fe] ratios correspond to larger 
Mgb indices. In the same figure we have plotted the glob- 
ular clusters data (Puzia et al. l2002p . the integrated bulge 



Table 3. The a/Fe bias in the Milky Way. 



z 


0.0004 


0.004 


0.008 


0.02 


0.05 


[a/Fe] 


0.22 


0.14 


0.1 


0. 


0. 



light (Puzia et al. i2002p , and the 65 elliptical galaxies of 
the sample of Papers I -I- II. The plot is analogous to that 
presented in Fig. 2 of TMB03 and Fig. 3 of K05. 

3.8. Calibration of an enhancement-independent index 

We have checked if the [MgFe] index, proposed by G93 
and revised by TMB03, is still a tracer of total metallic- 
ity independently of the a-enhancement. G93 defines the 
index 



[MgFe] = x/Mgh < Fe > 
with 

<Fe> = i(Fe5270-l-Fe5335), 



(23) 



(24) 



while TMB03 propose the new index [MgFe] , whose defi- 
nition is still given by Eq. ((23|) , although the < Fe > index 
is replaced by 



< Fe > 



(0.72 • Fe5270 + 0.28 • Fe5335). 



(25) 



In Fig. [6] we plot the indices Mgb, < Fe >, [MgFe], and 
[MgFe] as a function of [a/Fc] for our new SSP models. 
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Fig. 6. Lick indices as a function of [a/Fe] ratio 
at fixed solar metallicity. We plot the indices Mgb, 
< Fe >= i(Fe5270 + Fe5335), [MgFe], and [MgFe]' (see 
text for details) for a 12-Gyr-old SSP of our new a- 
enhanced models. 

In all the models, the Mgb index increases with the a- 
enhancement, while < Fe > decreases. On the other hand, 
the [MgFe] and [MgFe] indices remain essentially con- 
stant. More specifically, the fractional variations of [MgFe] 
and [MgFe] amount to ^ —0.004 and ~ 0.006, respec- 
tively, passing from [a/Fe]= to [a/Fe]= 0.5. 

3.9. Comparison with TMB03 SSPs 

We have compared our new SSP models with those of 
TMB03. Since both models are based on the fitting func- 
tions of Worthey et al. (J1994|) and Worthey & Ottaviani 
(|1997p . the possible differences between the models for 
solar-scaled compositions must be due to the different evo- 
lutionary tracks adopted. For a-enhanced compositions, 
an additional difference is the new index dependence on 
element abundance introduced in our models (see Sect. 
3.5). 

In Fig. [7] we plot the II/3 and [MgFe] indices pre- 
dicted by our models and those of TMB03 for solar total 
metallicity, ages in the range 1-15 Gyr, and [q;/Fc]=0, 0.5, 
plus the TMB03 models for [Z/H]=0.35 and [a/Fe]=0. In 
particular, we observe that our solar metallicity models 
are between the [Z/H]=0 and [Z/H]=0.35 TMB03 mod- 
els. Thus we expect our models to yield lower metallicities 
for the same line-strength index data. The comparison 
is also showed in Table 4. We compare our models with 
those of TMB03 and Tantalo & Chiosi (|2004a[ (TC04) for 
solar metallicity, [a/Fc]=0, and ages=5, 10 Gyr, for some 
commonly used indices. We observe that there is a sig- 
nificant difference in the metallic indices, in particular the 
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Table 4. Comparison with SSPs in the literature 





H/3 


Mg2 


Mgb 


Fe5270 


Fe5335 


[MgFe]' 


5 Gyr 


US 


1.95 


0.23 


3.60 


2.86 


2.59 


3.17 


TMB03 


1.99 


0.21 


3.18 


2.83 


2.52 


2.95 


TC04 


1.90 


0.23 


3.30 


2.90 


2.63 


3.06 


10 Gyr 


US 


1.62 


0.26 


4.08 


3.06 


2.79 


3.49 


TMB03 


1.70 


0.23 


3.58 


2.98 


2.66 


3.22 


TC04 


1.50 


0.27 


3.76 


3.12 


2.85 


3.38 



Mgb. Our models predict Mgb indices that are quite larger 
than TMB03, and this results in higher [MgFe]'. The TC04 
models are instead between our models and the TMB03 
ones. Bluer metallic indices are probably the results of 
hotter stellar tracks. In particular it is well known that 
the Girardi et al. (|2000p tracks used in the TC04 models 
predict warmer red giant branches (see, e.g., Bruzual & 
Chariot i2003j). 

To perform a more direct and quantitative comparison, 
we have considered the (H/3, Mgb, and <Fe>) indices for 
some selected TMB03 models of a given age, [Z/H], and 
[a/Fe], and have fitted those indices with our models. This 
allows us to directly compare the (age, [Z/H], and [a/Fe]) 
solutions provided by our models and those of TMB03 
for the same index triplet. Our test shows that there are 
significant differences in the derived stellar population pa- 
rameters. The major discrepancy concerns metallicity: our 
models provide metallicities that are systematically lower 
(up to 0.2 dex) than those predicted by the TMB03 mod- 
els. On the other hand, our derived ages tend to be older 
(10% higher at 10 Gyr). The [a/Fe] ratios are instead in 
good agreement. We only observe that for high [a/Fe] ra- 
tios 0.5), we tend to predict values that are larger by 
20 %. The same metallicity-age discrepancy has also been 
discussed by Denicolo' et al. (|2005bp when comparing the 
TMB03 models with those of Worthey (1994): fitting the 
line strength index data of the galaxy NGC 3379 with 
both the set of models, they derived a higher metallicity 
(A[Fe/i/] = 0.08) and a younger age (A age =-0.5 Gyr) 
with the TMB03 SSPs than with the Worthey ((TOM|) ones. 

4. Derivation of stellar population parameters 

4.1. The line-strength index diagnostic 

The use of absorption-line indices has proven to be a pow- 
erful tool in the derivation of ages and metallicities of un- 
resolved stellar populations (e.g., Worthey 1 1994p . To this 
purpose it is important to identify those absorption fea- 
tures, among the Lick indices, that possess a marked sen- 
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Fig. 7. Comparison of our SSP models with the models 
of TMB03 for solar total metallicity, ages in the range 1- 
15 Gyr, and [a/Fc]=0, 0.5. In the H/3 vs. [MgFe]' plane, 
our models are denoted with the solid line, the TMB03 
models with the short-dash line. The long-dash line is 
for the TMB03 models at [Z/H]=0.35. We observe that 
our solar metallicity models are between the [Z/H]=0 and 
[Z/H]=0.35 TMB03 models. 

sitivity either to abundance or age, that are not strongly 
affected by uncertainties in the measurement process, and, 
finally, that can be described by fairly robust models. 
Among the Balmer lines, which, measuring the presence of 
warm A-type stars, can be used as age indicators, the H/3 
index turns out relatively insensitive to abundance ratio 
variations, but it may be contaminated by nebular emis- 
sion. In reverse, higher order Balmer lines, less affected 
than Hz? by emission infilling (in particular the US less 
than the II7), are not independent from abundance ratio 
variations, as shown in Sect. 3. 

Line-strength indices, like the Mg and Fe indices, de- 
rived from metallic lines are mainly sensitive to element 
abundances, but they also depend on age. In particular, 
we have shown in Sect. 3.8 that the [MgFe] and [MgFe] 
indices are good tracers of total metallicity and are almost 
completely independent of a/Fc ratio variations. Thus, a 
three-dimensional space defined by a Balmer absorption 
line, the [MgFe] (or [MgFe] ), and a metallic index sensi- 
tive to a-enhancement (like Mgb) demonstrates itself to 
be a suitable diagnostic tool for the derivation of ages, 
metallicities, and [a/Fe] ratios. 

In Fig. m we plotted four projections of the spaces de- 
fined by (H/3, [MgFe]', Mgb), {R-jf, [MgFe]', Mgb), and 
(RSf, [MgFe]', Mgb). In the planes Mgb - [MgFe]', H/3 
- [MgFe]', BjF - [MgFe]', and B6f - [MgFe]', we super- 
posed the SSP models on the line-strength indices derived 
for our sample within re/8. Galaxies with different detec- 
tion levels of the [OIII]A5007 emission line are plotted us- 
ing different symbols. In this plot we show only 62 out of 
the initial sample of 65 galaxies. We have excluded: the 



Seyfert galaxy IC 5063 because of the presence of very 
strong emission lines, which prevented us from deriving 
reliable absorption line indices; NGC 6875, which lacks a 
measure of velocity dispersion (consequently the indices 
have not been corrected for this effect), and NGC 2962, 
which clearly shows calibration problems in the higher or- 
der Balmer indices and presents significant emission that 
can affect the H/3. These three galaxies will not be in- 
cluded in our results. In the Balmer indices vs. [MgFe] 
planes, models of constant age (almost horizontal lines) 
and of constant metallicity (almost vertical lines) are fairly 
well separated as a result of the markedly different sen- 
sitivity of the Balmer index to age and of the [MgFe] 
index to metallicity (see Worthev ll994p . The effect of a- 
enhancement is that of slightly strengthening the H/3 fea- 
ture. Though small, the variation at old ages corresponds 
to a few Giga years. On the contrary, [MgFe] is almost 
unaffected, due to the opposite response of Mg and Fe in- 
dices to a-enhancement. The H/3 vs. [MgFe] plane, where 
the enhanced SSPs are almost superposed on the standard 
ones, is in principle well suited to directly reading off ages 
and total metallicities almost independently of abundance 
ratio effects. In the H/3 vs. [MgFe] plane, but also in the 
other planes defined by high order Balmer lines, the data 
cluster in the lower right region of the diagram, between 
solar and twice solar metallicity models, and at old ages. 
A tail of objects extends at higher H/3 values and, as al- 
ready noticed by Longhetti et al. (|2000p . does not follow 
a line of constant metallicity. A large scatter in age is ob- 
served, and galaxies with equivalent-SSP ages of few Gyrs 
are present. 

There are, in particular, few galaxies that fall outside 
the models at very high H/3 values, even extrapolating the 
models at metallicities higher than Z— 0.05: NGC 1052, 
NGC 3136, NGC 5266, and NGC 2911. We have labeled 
these objects in all the four panels to better understand 
the origin of the discrepancy. In particular we want to 
know if the high H/3 values are caused by an overestima- 
tion of the emission correction or are rather due to in- 
trinsic properties. We observe that NGC 3136 has high 
values for all the three Balmer indices. We note in partic- 
ular that while the H/3 has been corrected for emission, 
H'jF and HSf have not been. Thus the plotted H7j;^ and 
HSf are lower limits, and by applying the (small) emis- 
sion correction to them, we would get even higher values. 
Given the high quality and the high S /N of the spectrum 
also, we conclude that the observed strong Balmer lines 
for NGC 3136 are a real effect. As noted by Bressan et 
al. (1996) and Longhetti et al. (2000), a galaxy in which 
a recent burst of star formation is superimposed on an 
old stellar population can fall above the H/3 - [MgFe] re- 
gion occupied by the SSPs. The reason is that if in the 
secondary burst the gas metallicity was higher than the 
metallicity of the old population, the galaxy moves verti- 
cally toward the top of the H/3 - [MgFe] plane. Thus, it is 
likely that the position of NGC 3136 in the H/3 - [MgFe]' 
plane is due to the occurrence of secondary bursts of star 
formation. 
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Fig. 8. Lick line-strength indices (re/8 aperture) of our sample compared with SSP models with metallicities Z= 
0.008,0.02, and 0.05, ages in the range (2-15) Gyr, and [a/Fe]=0 (solid line) and 0.4 (dotted line). We have plotted 
H/3 - [MgFe]' (top left panel), Mgb - [MgFe]' (top right), H7F - [MgFe]' (bottom left), and RSp - [MgFe]' (bottom 
right) planes. The metallicity increases for increasing [MgFe] values, while the age increases for decreasing Balmer 
index. Vertical lines in the Balmer index vs. [MgFe] plane are models of constant metallicity, while horizontal lines 
are models of constant age. The data have been plotted with different symbols according to their emission properties: 
open circles for [OIIIA5007] emission detected under la, triangles for emission between 1 and 2 a, and squares > 2a. 
The galaxies IC 5063, NGC 6875, and NGC 2962 have been excluded from the original sample of 65 galaxies. 



The situation seems to be different for NGC 1052, 
NGC 5266, and NGC 2911. They all show high H/3 val- 
ues, but H7F and HSp fall within the models (NGC 5266 
and NGC 2911) or are too low (H^f for NGC 1052). The 
observed inconsistency among the Balmer lines seems to 
be due to an overestimation of the H/3 emission, with con- 
sequent overcorrection of the index. Given the accurate 
subtraction of the continuum around Ha performed to de- 



termine the emission (see Papers I and II), the only viable 
explanation for the high 11(3 emission is the presence of sig- 
nificant dust extinction: in this case the FHa/FII/3 ratio 
is higher than 2.86, which is the value adopted through- 
out our analysis. Indeed HST observations have revealed 
dust absorption features in the center of NGC 1052 (van 
Dokkum & Franx[inMl Gabel et al. lMID)) : NGC 2911 is 
a disk dominated SO with an important dust component 
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(Michard & Marchal [HMD , and NGC 5266 is an Hl-rich 
elliptical galaxy with a dust ring along the minor axis 
(Varnas et al. [TM71 Morganti et al. [TM7|) . 

The index errors displayed in Fig. 8 have been ob- 
tained adding a Poissonian fluctuation to each spectrum 
and generating a set of 1000 Monte Carlo modifications 
(see Papers I and II). The error in the H/3 index does not 
include the systematic uncertainties in the emission cor- 
rection due to the choice of a particular template for the 
continuum subtraction and of a FHa/FH/? ratio. Thus, we 
warn the reader that the H/3 error bars plotted in Fig. 8 
may be underestimated. 

In the Mgb - [MgFe] plane, models of constant age and 
metallicity are almost degenerate, while the effect of the 
enhancement is well separated. We observe in particular 
that the bulk of data is confined between [a/Fe]= and 
[a/Fe]= 0.4. Furthermore, there is a trend for galaxies 
with larger [MgFe] values to lie closer to the enhanced 
models than to the standard ones, suggesting a correlation 
between metallicity and a/Fe enhancement. The three- 
dimensional spaces described in this section will be used 
in our following analysis to derive ages, metallicities, and 
[a/Fe] ratios for our sample of galaxies. For this purpose 
we have devised a simple but robust algorithm, which we 
describe in the next section. 

4.2. The probability distribution of the physical 
parameters 

A SSP model of given age, metallicity, and a/Fe ratio cor- 
responds to a point in a multi-dimensional space defined 
by the line-strength indices. If we exclude ages older than 
~ 17 Gyr, where a blueing of the indices may be intro- 
duced by the presence of Horizontal Branch (HB) stars 
even at high metallicity, there is a unique correspondence 
between the index multiplet and a point in the age, metal- 
licity, and [a/Fe] parameter space. In principle, one could 
chose three arbitrary indices and, by comparison with SSP 
models, derive the (age, Z, [a/Fe]) solution. Nevertheless, 
as discussed in the previous section, some combinations of 
indices among the Lick set are more useful than others to 
disentangle the age-metallicity degeneracy, and thus turn 
out to be particularly suitable for deriving the stellar pop- 
ulation parameters. More specifically, for our analysis we 
select the three dimensional spaces defined by (II/3, <Fe>, 
Mgb), (H7F, <Fe>, Mgb), and (RSf, <Fe>, Mgb). Once 
the index 3D space has been defined, the stellar popu- 
lation parameters are derived according to the following 
procedure. 

Let us call (xq, j/oj -^o) the observed point in the in- 
dex space and {a.j;„ , <Jy„ , a^a ) the associated error. For each 
point (t, Z, a) in the (age, metallicity, a/Fe enhancement) 
space (G space), we can define the probability density that 
the generic SSP model ix(^t,z,a),y{t,z,a), Z(t.z,a)) is the so- 
lution to the observed line-strength indices. In particular, 
if the errors are Gaussian, the probability density assumes 
the expression: 
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Equation (|26p allows us to construct a probability density 
map Pt,z,a defined on the three dimensional grid {t, Z, a). 

In Fig.[9l we show the contour levels in the (t, Z) plane, 
enclosing the solutions that are within one a (in the three 
indices) from the observed values for some galaxies of our 
sample, as an example. Each panel refers to a different 
galaxy, and the three curves refer to the solutions obtained 
with the different triplets adopted, i.e., II/3, <Fe>, Mgb; 
H7i?, <Fe>, Mgb; and YiSp, <Fe>, Mgb. The most prob- 
able solution in each index triplet is obtained by finding 
the maximum of the Pt,z,a function. 

Figure [9] shows that measurement errors combine with 
the relative index sensitivity to age and metallicity vari- 
ations (the age metallicity-degeneracy). Levels of con- 
stant probability define a narrow elongated region of anti- 
correlation between age and metallicity solutions. Some 
authors have shown this effect by using Montecarlo simula- 
tions to generate a mock catalog of synthetic galaxies (e.g., 
Thomas et al. I2005P . Notice that the age-metallicity de- 
generacy introduces an asymmetry in the errors of the so- 
lution, as a Gaussian error in the measured indices trans- 
lates into a non-Gaussian error in the solution. We have 
used the three triplets to minimize spurious observational 
effects in the age indicator H/3, which may be affected by 
emission filling. For each galaxy we combined only those 
index triplets that provide more reliable results, depend- 
ing on the strength of the present emission (see below). 

For most of the sources, the contour levels of the three 
index spaces possess a clear region of intersection. This 
indicates that models predict consistent results. However, 
there are galaxies for which the three G spaces do not 
have common solutions (within one-cr, if we consider two- 
(7, then there are always common solutions, but the uncer- 
tainty is large). In particular, for some galaxies, H/3 pro- 
vides solutions that are really detached from those pro- 
vided by the other two age indicators. As discussed in 
Sect. 4.1, discrepancies between the H/3 and the higher 
order Balmer lines could be due to problems in the deriva- 
tion of the emission correction. 

For each 3D index space, the solution is computed 
adopting a subspace (G*) defined by Pt,z,a > / x Pmax, 
and performing an average on the (G*) values weighted 
for the probability density: 



" Jj/g. P{t, Z, a) dt dZ da 



JJJ^, Z P{t, Z, a) dt dZ da 
JJJ P(t, Z, a) dt dZ da 



(27) 
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Fig. 9. Example of contour levels enclosing solutions that are within one a from the observed value. Black solid lines 
are for (H/3, <Fe>, Mgb) triplets, red dotted lines for the (H7i?, <Fe>, Mgb) triplets, and blue dashed lines for the 
(H(5i?, <Fe>, Mgb) triplets. Diamonds, triangles, and squares refer to the most probable solution for each triplet, while 
crosses indicate the average of the logarithmic values of the solutions. 



///^,aP(t,Z, 



dt dZ da 



JJ/q. P{t, Z, a) dt dZ da 



More specifically, we consider the subspace of solutions 
obtained with f=0.9. The uncertainties on the solutions 
are computed considering the whole G space: 



ffj^it - tf,y P{t, Z, a) dt dZ da 



IIIg ^' '^^ 



IIlGiZ 



Z^,f P{t,Z,a) dt dZ da 



IIIg ^(^' ^' ") '^^ 
Illaia - a^)2 P{t, Z, a) dt dZ da 
Jj/p P{t, Z, a) dt dZ da ' 



(28) 



or absent emission {\EWem{H(3)\ < 0.2), there is a good 
agreement among the stellar population parameters de- 
rived with the use of the three different Balmer indices. 
For larger emissions, the results obtained with the H/3 in- 
dex deviate from the H7 and predictions however, with 
a tendency toward younger ages. Finally, for strong emis- 
sions {\EWem^{HP)\ > 0.6), the infilling of the H7 starts 
to become significant and the index, which has not been 
corrected for this effect, systematically provides very old 
ages. In this last case, the US index, that presents the low- 
est emission contamination, is perhaps the most reliable 
one. Thus we combine the information from the different 
index triplets, performing a weighted mean over those N 
indices which, according to the galaxy emission properties, 
provide a reliable solution as explained above: 



The analysis, repeated for the triplets (H/?, <Fe>, Mgb), 
(H7F, <Fe>, Mgb), and {HSp, <Fe>, Mgb), provides 
three different partial solutions, tf^^i, and q;^,^ for i=l, 
3, with the corresponding uncertainties. 

At this point we compare the results provided by the 
different index triplets and check for consistency within 
the errors. We observe that for those galaxies with low 
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[a/Fe] (Mgb) 

Fig. 10. Comparison of age, metallicity, and a- 
enhancement values derived using Mg2 instead of Mgb 
index. 



where the weights are derived from the uncertainties on 
the parameters: 

Wt,, = -^■,Wz,^ = ^;W^,^ = ^. (30) 

'^Z,^ 



5. Stellar populations in early-type galaxies 

Through the procedure described in Sect. 4.2, we have de- 
rived stellar population parameters (age, metallicity, and 
[a/Fe] ratio) for the apertures and gradients of the 62 
galaxies plotted in Fig. 8. At the base of our analysis is 
the assumption that early-type galaxies contain almost co- 
eval stellar populations and may be well described by only 
one SSP. This is the most direct method to derive hints 
about the correlations between stellar population param- 
eters and galaxy properties, like mass and environment. 
In the following analysis, we will focus both on the cor- 
relation between central galaxy populations and velocity 
dispersion/environment and on the analysis of the radial 
gradients. 

5.1. Ages, metallicities, and a/Fe in the central region 

We provide ages, metallicities, and [a/Fe] ratios derived 
from the three G-spaces (H^, <Fe>, Mgb), (H7F, <Fe>, 
Mgb), and {HSp, <Fe>, Mgb) for the apertures of the 
sample in Table 5 to 7. Column (1) gives the galaxy iden- 
tification name. Col. (2) the aperture number (see Papers I 
and II for the aperture radii), and Cols. (3), (4), and (5) 
provide the age, the metallicity mass fraction Z {Z = 0.018 
is the solar value) and the [a/Fe] ratio respectively. The 
complete tables are given in electronic form. In Table 8 
we provide the final combined solutions for a re/8 aper- 
ture, obtained through Eqs. (26) and (27). Columns (1) 
and (2) give the galaxy identification name and the central 
velocity dispersion respectively. Columns (3), (4), and (5) 
provide the age, the metallicity mass fraction Z (Z = 0.018 
is the solar value), and the [a/Fe] ratio, respectively. 

Before continuing the discussion, we remark that our 
results are fairly independent of the choice of cither mag- 
nesium index. Figure [10] compares the results obtained 
using the Mg2 index instead of the Mgb index. We no- 
tice a very good consistency in the derived ages, while 
the metallicity and a-enhancement would be only slightly 
higher adopting Mg2 instead of Mgb. However, the size of 
the latter differences will not affect our conclusions. 

The distributions of age, metallicity, and enhancement 
within re/8 are shown in Figure 1111 In each panel, we 
present the results obtained from the three different G 
spaces. We plot both the distributions for the total sample 
and the partial distributions for lenticular and elliptical 
galaxies. For all three index triplets, the age distribution 
shows a large spread (1 Gyr< age< 15 Gyr). From the 
H/3 and H7 indices, a bimodality, with a group of old ob- 
jects (ages > 7 Gyr), mainly composed by ellipticals, and 
a tail of young objects, mainly lenticulars, appears. We 
also observe that H/3 tends to provide a larger number of 
very young galaxies than II7 and US . This effect is likely 
the result of the high uncertainty in the applied emission 
correction. However, we stress again that the method we 
adopted to combine the results from the different indices 
(see Sect. 4.2) allows us to clean our sample from spuri- 
ous effect and obtain unbiased results for the final ages. 
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metallicities, and [a/Fe] ratios. The average ages for the 
whole sample, E and SO obtained by properly combining 
the information from the three index triplets, are 8, 8.7, 
and 6.3 Gyr, respectively. 

The metallicity distribution shows a broad peak, lo- 
cated between 0<[Z/H]<0.3. There is a good agreement 
among the values derived from the analysis of the differ- 
ent index triplets, with the H(5 distribution being peaked 
at slightly lower metallicities, which is an evident effect of 
the age-metallicity degeneracy. The average metallicities 
for the whole sample, E and SO, are 0.21, 0.22, and 0.19, 
respectively. Finally, the [a/Fe] ratio presents a maximum 
at ~ 0.22 with a narrower distribution than those of age 
and metallicity. The average ratios for the whole sample, 
E and SO, are 0.21, 0.23, and 0.17, respectively. 

5.1.1. Correlation with central velocity dispersion 

Figure [12] shows selected Lick indices measured in the cen- 
tral aperture (r < re/8) for the 65 galaxies of our sample 
as a function of central velocity dispersion. We observe 
that the separation between the two classes (those galax- 
ies classified as E according to RC3 and those classified as 
SO) is quite well marked in velocity dispersion, with very 
few E galaxies populating the region below CTc ~ 180 km/s. 
The fit to the total sample with coefficients is given in the 
upper label. We also show the percentage index variation 
along the regression when the velocity dispersion changes 
from CTc = 300 km/s to CTc = 100 km/s. Separate fits are 
made for the SO and E galaxies. 

It is clear that metallic indices show a well-established 
positive trend with Uc (with a significantly shallower vari- 
ation for the Fe indices). On the contrary, Balmer indices 
decrease with velocity dispersion. The observed behavior 
suggests metallicities or ages increasing with galaxy mass, 
but for the quantitative results we refer to the following 
comparison with SSP models. 

To seek for possible correlation of the derived stel- 
lar population parameters with central velocity disper- 
sion/galaxy mass, we plot the final solutions derived 
within a re/8 aperture (see Table 8) as function of ctc 
in Fig. 1131 In the top, central, and bottom panels, we 
plot ages, metallicities, and [a/Fe] ratios, respectively. 
Figure [13] shows a trend of increasing age with central 
velocity dispersion. Considering the separate fits for el- 
lipticals and Icnticulars, we obtain that the trend is pos- 
itive for lenticulars and negative for ellipticals. However, 
a Spearman rank order test provides no proof for the ex- 
istence of a correlation between age and CTc — for both 
the whole sample and the E and SO subsamples. The sep- 
arate fits suggest that on average lenticulars are younger 
than ellipticals. The former, however, shows a large scat- 
ter. There is also a group of young E galaxies that deviate 
significantly from the main trend, at CTc ~250 km s~^. 
The middle panel shows that there is a significant posi- 
tive trend of metallicity with Cc, with no significant dif- 
ference in the slope between the two morphological types. 



However, lenticular galaxies appear on average more metal 
rich than ellipticals at the same velocity dispersion. A 
Spearman test provides p-levels of significancqj of 0.0001, 
0.0002, and 0.07 for the whole sample, E and SO, re- 
spectively, confirming that there is strong correlation be- 
tween the parameters. The same result is found for the 
a-enhancement, which again increases with central veloc- 
ity dispersion. For the whole sample, the Spearman test 
results in a p level of 0.0008. The absence of a clear trend 
of age with ctc, and the presence of strong correlations for 
metallicities and a-enhancement suggest that the galaxy 
gravitational potential must mainly affect the chemical 
enrichment process of the galaxy. 

5.1.2. Correlation with environment 

A measure of the richness of the environment, Pxyz, sur- 
rounding each galaxy in galaxies Mpc~^, is reported by 
TuUy (1988) for 47 galaxies (73% of the whole sample) 
(see Table 1). Even though our sample is prevalently made 
of field galaxies and we are far from the densities found in 
rich galaxy clusters, we may investigate how (eventually) 
the environment affects the trend of the physical param- 
eters with velocity dispersion. In Fig. [TJ] we plot the de- 
rived age, metallicity, and enhancement values within an 
re/8 aperture as a function of the richness parameter Pxyz- 
Analogously to Fig.[T3] E and SO galaxies are denoted, and 
a linear fit is performed for all the galaxies and for the E 
and SO subsamples. The upper panel indicates a positive 
trend of age with the richness parameter (the Spearman 
test provides a 85 % probability for a correlation to ex- 
ist). As for what concerns metallicity and a-enhancement, 
there are no clear relations with p. The fitted relations to 
the data suggest that galaxies in lower density environ- 
ments are slightly more metal rich and less enhanced in a 
-elements than galaxies in richer environments. However, 
a Spearman test does not provide any evidence for the 
existence of significant correlations. 

The relation of age with environment is particularly 
interesting, in view of the fact that we have not found 
any trend with velocity dispersion. It is surprising that, in 
spite of containing more than 40% of the sample with de- 
termined Pxyz, the region above \og{pxyz) >-0.4 does not 
contain galaxies younger than 4 Gyr. If we consider only 
galaxies with ages older than 4 Gyr, the resulting fit is 
fiat, indicating that the age-environment relation is actu- 
ally due to the presence of very young objects in the poorer 
environments. In fact, in such environments we have the 
contemporary presence of objects with both high and low 
velocity dispersions. It is also clear, after comparison with 
Fig. [131 that the clump of young E galaxies at Cc —250 
km s~^ is due to objects located in a very low density 
environment. 

Figure [14] also suggests that the relation with velocity 
dispersion may be " disturbed" by the presence of such ob- 

^ The p level of significance represents the probability that 
the relation found between the variables is not true. 
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Table 8. Average ages, metallicities, and a/Fe ratios (re/8 apert 



Ident. 


km 


Age 
Gyr 




Z 






:/Fe] 


NGC 128 


183 


9.7 ± 1.7 


0.024 


± 


0.004 


0.16 


± 


0.03 


NGC 777 


317 


5.4 ± 2.1 


0.045 


± 


0.020 


0.28 


± 


0.10 


NGC 1052 


215 


14.5 ± 4.2 


0.032 


± 


0.007 


0.34 


± 


0.05 


NGC 1209 


240 


4.8 ± 0.9 


0.051 


± 


0.012 


0.14 


± 


0.02 


NGC 1297 


115 


15.5 ± 1.2 


0.012 


± 


0.001 


0.29 




0.04 


NGC 1366 


120 


5.9 ± 1. 


0.024 


± 


0.004 


0.08 




0.03 


NGC 1380 


240 


4.4 ± 0.7 


0.038 


± 


0.006 


0.24 


± 


0.02 


NGC 1389 


139 


4.5 ± 0.6 


0.032 


± 


0.005 


0.08 


± 


0.02 


NGC 1407 


286 


8.8 ± 1.5 


0.033 


± 


0.005 


0.32 


± 


0.03 


NGC 1426 


162 


9.0 ± 2.5 


0.024 


± 


0.005 


0.07 


± 


0.05 


NGC 1453 


289 


9.4 ± 2.1 


0.033 


± 


0.007 


0.22 




0.03 


NGC 1521 


235 


3.2 ± 0.4 


0.037 


± 


0.006 


0.09 


± 


0.02 


NGC 1533 


174 


11.9 ± 6.9 


0.023 


± 


0.020 


0.21 


± 


0.10 


NGC 1553 


180 


4.8 ± 0.7 


0.031 


± 


0.004 


0.10 


± 


0.02 


NGC 1947 


142 


5.9 ± 0.8 


0.023 


± 


0.003 


0.05 




0.02 


NGC 2749 


248 


10.8 ± 2.3 


0.027 


± 


0.006 


0.25 


± 


0.04 


NGC 2911 


235 


5.7 ± 2.0 


0.034 


± 


0.019 


0.25 


± 


0.10 


NGC 2974 


220 


13.9 ± 3.6 


0.021 


± 


0.005 


0.23 


± 


0.06 


NGC 3136 


230 


1.5 ± 0.1 


0.089 


± 


0.004 


0.36 


± 


0.02 


NGC 3258 


271 


4.5 ± 0.8 


0.047 


± 


0.013 


0.21 




0.03 


NGC 3268 


227 


9.8 ± 1.7 


0.023 


± 


0.004 


0.34 




0.04 


NGC 3489 


129 


1.7 ± 0.1 


0.034 


± 


0.004 


0.05 


± 


0.02 
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5.8 ± 0.8 


0.034 


± 


0.004 


0.17 


± 


0.02 


NGC 3607 


220 


3.1 ± 0.5 


0.047 


± 


0.012 


0.24 




0.03 


NGC 3818 


191 


8.8 ± 1.2 


0.024 


± 


0.003 


0.25 


± 


0.03 


NGC 3962 


225 


10.0 ± 1.2 


0.024 


± 


0.003 


0.22 


± 


0.03 


NGC 4374 


282 
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± 
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± 


0.002 
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± 
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± 


0.05 
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± 
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0.19 


± 


0.03 


NGC 6876 
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± 


0.003 
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± 
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± 


0.006 
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± 
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0.010 


0.15 
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Table 8 Continue. Average ages, metallicities, and a/Fe ratios (re/8 aperture) 



Ident. 




Age 


Z 


[a/Fe] 




km s^^ 


Gyr 






IC1459 


311 


8.0 ± 2.2 


0.042 ± 0.009 


0.25 ± 0.04 


IC2006 


122 


8.1 ± 0.9 


0.026 ± 0.003 


0.12 ± 0.02 


IC3370 


202 


5.6 ± 0.9 


0.022 ± 0.004 


0.17 ± 0.04 


IC4296 


340 


5.2 ± 1.0 


0.044 ± 0.008 


0.25 ± 0.02 



Notes: the values are obtained combining the H/3, H7, H5, Mgb, and <Fe> indices as described in Sect. 4.2. 




Age (Gyr) [Z/H] [a/FeJ 



Fig. 11. Distribution of ages (left panel), metallicities (central), and [a/Fe] ratios (right) measured at re/8, derived 
from the (H/3, <Fe>, Mgb), (H7F, <Fe>, Mgb), and (RSf, <Fe>, Mgb) index spaces. In the plot we have included 
only the galaxies for which the solution is derived within 1 a from the observed indices. The solid line refers to the 
total sample, while the dotted and dashed lines are for SOs and Es, respectively. For each distribution, the used index 
triplet is provided in the label inside the box. 



jects populating very low density environments. In Fig. 1151 
we have redrawn the fit with velocity dispersion separat- 
ing objects according to their richness parameter, instead Log(age/Gyr) = -0.15 + 0.43 x log((Jc) (Total : r = 0.2) 
of morphological type. A linear fit is performed for all the Log(age/Gyr) — —0.51 + 0.56 x log((jc) (LDE : r = 0.24) 
galaxies, while fits to the galaxies in high and low den- Log(age/Gyr) = 0.86 + 0.03 x log(o-c) (HDE : r = 0.01) 



sity environments are also performed separately. We notice 
that at values log{pj:yz) ~ -0.4, we may identify small asso- 
ciations of galaxies, such as pairs and poor galaxy groups, 
in the TuUy (|1988|) catalogue. On the upper side, above 

log(p^j,^) « 0.0 typical ofthe outskirts ofVirgo cluster, we Log(Z/ZQ) = -0.99 4- 0.48 x log(ac) (Total : r = 0.34) 
count only four objects, and we stress again that our rich Log(Z/ZQ) = — 0.85 -I- 0.44 x log((Jc) (LDE : r = 0.3) 
environment is not representative of rich clusters. Again Log(Z/ZQ) = -1.36 -I- 0.64 x log((Jc) (HDE : r = 0.47) 
there is no clear trend of age with velocity dispersion, but 
now galaxies in relatively richer environments are more 
tightly clustered around the average relation. The follow- 
ing correlations are derived for age, metallicity, and a/Fe 

enhancement, respectively: [a/Fe] = -0.74 -I- 0.41 x log((Jc) (Total : r = 0.53) 

[a/Fe] = -0.79 + 0.43 x log((Jc) (LDE : r = 0.52) 



(31) 



(32) 
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Fig. 12. Selected Lick indices, measured at r < re/8, for the total sample of 65 galaxies as a function of log(crc/300), 
where CTc is the central velocity dispersion. Triangles and diamonds indicate E and SO galaxies, respectively. The dashed, 
dotted, and solid lines mark the linear fit obtained for SO galaxies, E galaxies, and the total sample, respectively. For 
each index, the linear fit to the total sample is labeled above each panel. 



[a/Fe] = -0.87 + 0.47 x log(crc) (HDE : r = 0.58), 

(33) 

where LDE and HDE indicate galaxies populating low and 
high density regions, respectively, and r is the linear cor- 
relation coefficient. 

At 200 km s^i, the average age of galaxies inhabit- 
ing HDE is of about 8.5 Gyr, while those in LDE is ~6 
Gyr. As far as the relation between metallicity and CTc is 
concerned, we notice that there is almost no dependence 
on the environment. The different relations are almost su- 
perimposed, with the only noticeable difference being that 
the galaxies in HDE are less dispersed than those in LDE. 
The same is valid for the enhancement. 

Summarizing, combining the information of Fig. [15] 
with the behavior of the stellar population parameters 
with environment, we argue that the metal enrichment is 
essentially dominated by the gravitational binding of the 
parent galaxy. In reverse, the way in which the popula- 
tions are assembled (the growth of baryons within galax- 
ies) seems modulated by the environment, with galaxies in 
richer environments being on average older than galaxies 



in less rich environments. Galaxies in more rich environ- 
ments show a lack of very young members. Since environ- 
ment seems to have no effect on the metallicity-cTc and 
[a/Fe\~(Tc relations, we argue that very young members in 
the lowest density regions are more likely due to rejuve- 
nation episodes than to more prolonged star formations. 

5.2. Stellar population gradients 

Gradients within galaxies can provide further insight on 
the process of formation and evolution of early-type galax- 
ies. Our sample shows the presence of clear gradients in 
the observed line-strength indices. Gradients have been 
extracted from a linear fit to the index values measured in 
the four radial regions < ?' < re/16, re/16 < r < re/8, 
re/8 < r < re/4, and re/4 < r < re/2. 

In Fig. [Tni we show as an example the observed in- 
dex gradients as a function of CTc for some selected in- 
dices. Analogously to Fig. [131 we have used different sym- 
bols for SO and E galaxies, and different lines to denote 
the fits to SO and E separately, and the fit to the to- 
tal sample. The figure, besides showing the presence of 
eventual gradients of indices, highlights whether a corre- 
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Fig. 13. Ages, metallicities, and [a/Fc] ratios, measured at re/8, vs. the central velocity dispersion, logdc- The plotted 
stellar population parameters are derived from the combined analysis of the (H/3, <Fe>, Mgb), (ajF, <Fe>, Mgb), 
and (H(5f, <Fe>, Mgb) index spaces (see text for details). Triangles and diamonds denote, respectively, E and SO 
galaxies. The solid line is the linear fit performed to all the galaxies, while dashed and dotted lines are the best fit to 
Es and SOs subsamples. 



O 
ISI 

tS] 



1.2 

1.0 




■AA. 

^- 

A 




A A 




A 






-ffl- 




0.8 


f ■ -a 


— !r-^ 


^^^^ 














0.6 
0.4 
0.2 


TA 






O 

o 

o 














0.8 
0.6 


-A 




















0.4 
0.2 


— a 
: A 

^- - «• - - 




A 

— 


CO 


A 












0.0 
0.2 


r o 


O 




o 

A 


A A^ 








s 




0.3 


:A 


A 
^- 


A 
ZA 


o 


A 

A-sO 


A 




A 






0.2 

0.1 




= -A^^ 

r * o 


A 


.A. .4 

A A 
O 


o 

o 


A 

O 

A 




O 








-1 


.0 






0.5 


0.0 






0.5 




1 



Log(p) 



Fig. 14. Ages, metallicities, and [a/Fe] ratios, measured at re/8, vs. the density of the environment, log(p), in galax- 
ies Mpc""^ (TuUv [T9881 see Table 1). Triangles and diamonds denote, respectively, E and SO galaxies. The solid line 
is the linear fit performed to all the galaxies. Dashed and dotted lines are the best fit to E and SO subsamples. The 
dot-dashed line is the fit for the galaxies older than 4 Gyr. 
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Fig. 15. Ages, metallicities, and [a/Fc] ratios, measured at r-g/S, vs. the central velocity dispersion, logCc. In this 
figure triangles denote galaxies with a relatively higher richness parameter, log(p)>- 0.4. Diamonds denote the other 
objects, while crosses denote the total sample. The solid line is the linear fit performed to all the galaxies. The dashed 
and the dotted lines refer to galaxies in high and low density environments, respectively. 



lation between the gradient and the gravitational poten- 
tial of the galaxy exists. From Fig. [111 we notice that in 
general metallic indices decrease from the center of the 
galaxy outwards, while Balmer indices (i.e., H/3, USa) in- 
crease. Such a behavior suggests that either metallicity 
or age decreases from the central regions of the galaxy 
to the periphery. As far as the modulation of the gradi- 
ents with velocity dispersion is concerned, we notice that 
larger gradients tend to correspond to lower CTc values 
(at least for the metallic indices), suggesting that stel- 
lar population gradients tend to flatten with increasing 
galaxy mass. However, these are only qualitative con- 
clusions, and the quantitative derivation of the stellar 
population gradients requires the comparison of the ob- 
served indices with SSP models. Age, metallicity, and 
[a/Fe] values derived for the gradients (0 < r < re/16 
,re/16 < r < re/8, r^/S < r < re/4, re/4 < r < re/2) 
are provided for the three G-spaces (H/3, <Fe>, Mgb), 
(H7F, <Fe>, Mgb), and (RSp, <Fe>, Mgb) and for the 
final solution (Eq. [29]) in Tables 9 to 12. The complete 
tables are given in electronic form. For a more detailed 
presentation of the results with complete figures for the 
stellar population gradients, we refer to Annibah (|2005p . 
Gradients are expressed as 51og(Age/Gyr)/(51og(r/re), 
(51og(Z)/(51og(r/re), and (5[Q!/Fe]/(51og(r/re), respectively. 

We plot the distributions of age, metallicity, and a- 
enhancement gradients in Fig. [T7| . The age and [a/Fe] 
gradients present a large dispersion. The age gradients are 



symmetrically distributed around zero, while the [a/Fe] 
gradients show a low peak at negative values. The metal- 
licity gradient distribution is the only one that presents a 
clear peak at negative values. We derive an average metal- 
licity gradient of -0.21. 

As was done for the stellar population parameters de- 
rived within a re/8 aperture, we investigate possible trends 
of the gradients with the richness parameter p (Fig. [T8)) 
and the central velocity dispersion fFig. [T9|) . In Fig. [TSlE 
and SO galaxies are denoted. A linear fit is performed for 
all the galaxies, and fits are also performed for the E and 
SO subsamples. Like in Figure [T31 a fit to all the galax- 
ies older than 4 Gyr is performed. We notice that for the 
gradients there is not the dichotomy, revealed by Fig. [T4l 
between objects of relatively high and relatively low rich- 
ness parameters. However, there are hints of the existence 
of trends in the gradients with the richness parameter. The 
most evident trends are revealed by the chemical enrich- 
ment path. The metallicity gradient is at values around - 
0.25 (in good agreement with, e.g., Davies et al. ll993p . but 
seems to vanish at increasing richness. The tendency seems 
even more clear if only ellipticals are considered. The gra- 
dient of a-enhancement is slightly negative at very low 
density and increases with the richness parameter. The 
dispersion is, however, larger than in the case of the metal- 
licity gradient, and, in spite of the trend suggested by the 
linear fits, according to a Spearman rank order test, there 
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Fig. 16. Gradients of Lick indices computed as 5I/Slog{r) for the total sample of 65 galaxies as a function of 
log(crc/300), where CTc is the central velocity dispersion. Triangles and diamonds denote, respectively, E and SO galax- 
ies. The dashed, dotted, and solid lines mark the linear fit obtained for SO galaxies, E galaxies, and the total sample, 
respectively. For each index, the linear relation that fits the total sample is labeled above each panel. 



is no significant correlation between the gradients and the 
density parameter p. 

In Fig. [19] we consider the relation between gradients 
and central velocity dispersion once, as in Fig. 1151 we have 
separated galaxies according to their richness parameter. 
There is no significant trend of the age gradient with ve- 
locity dispersion for the whole sample. Again, the most 
important trends are those of metallicity and a/Fe- en- 
hancement. For the whole sample, the metallicity gradient 
shows an increasing trend with velocity dispersion (larger 
negative gradients for less massive systems); a Spearman 
test provides a significance level of 0.07 (93% probability 
that a correlation exist). The correlation is very strong 
if only E galaxies are considered {p = 0.009). As for 
the [a/Fe] gradient, it is slightly positive and constant 
for HDE galaxies ((5[a/Fe]/(51og(r/re)-0.05). LDE galax- 
ies have preferentially negative gradients, with all galaxies 
above 200 km/s lying below the relation defined by HDE 
galaxies. There are essentially no differences if we consid- 
ered the sample of E and SO galaxies separately (the plot 
is not shown here). 



6. Discussion 



How and when early-type galaxies formed is still a matter 
of debate. Two alternative formation scenarios have been 
proposed in the past. In the monolithic one, the spheroidal 
component forms by the gravitational collapse of a gas 
cloud with considerable energy dissipation (Larson [l974bl 
Arimoto & Yoshii 11987^ . As a result of the rapidity of 
this collapse, the bulk of stars in ellipticals should be rel- 
atively old. Unfortunately, this scenario is not supported 
by detailed hydrodynamical calculations: any angular mo- 
mentum initially present in the overdensity region pre- 
vents the rapid collapse of the gas. Small disky entities 
are formed first, and only after subsequent major merging 
events may they give rise to the spheroidal component of 
the galaxies (Toomre ll9771 Kauffmann et al. 119931 Baugh, 
Cole & Frenk 19961 Steinmetz & Navarro 12002]) . The star 
formation process is delayed in this so-called hierarchical 
scenario, which then predicts the presence of significant 
intermediate-age stellar populations. 
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Fig. 17. Distributions of the age, metallicity, and [a/Fe] gradients computed as Slog{Age/Gyr)/6log{r/re), 
^ log(Z)/(51og(r/re), and (5[a/Fe]/i51og(r/re), respectively. The dotted vertical line separates negative from positive 
gradient values. 



XJ 



1.0 - 

0.5 ~- 

0.0 M 

-0.5 ~- 

-1.0 ~- 
0.5 ~ 

0.0 - 

-0.5 - 

-1.0 - 
1.0 ~ 

0.5 r 
0.0 |=-| 

-0.5 r 

-1.0 L_ 
-1.0 



o 

A 

A AA 



A O 



-0.5 



A A 



A- A ^ • 



A <P 



0.0 

Log(p) 



0.5 



1.0 



Fig. 18. Age, metallicity, and [a/Fe] gradients against the density parameter, log(p). Triangles and diamonds represent, 
respectively, E and SO galaxies. The solid line is the linear best fit performed to all the galaxies. Dashed and dotted 
lines are the best fit to the E and SO subsamples. The dot-dashed line best fits all galaxies older than 4 Gyr. 



The study of absorption line indices in local early-type 
galaxies has proven to be one of the most powerful diag- 
nostics to constrain star formation history and to trace 
star evolution over time. Besides "direct" age estimates, 
indirect evidences of the duration of the star formation 
process can be analyzed through the study of the chemi- 
cal enrichment pattern, and more specifically of the a/Fe 
enhancement. While the a-elements O, Na, Mg, Si, S, Ar, 
Ca, and Ti (particles that are built up with a-particle 
nuclei) plus the elements N and Ne are derived mainly 
by Type II supernova explosions of massive progenitor 
stars, a substantial fraction of the Fe-peak elements Fe and 



Cr comes from the delayed exploding Type la supernovae 
(Nomoto et al. lT984l Woosley & Weaver flMSl Thielemann 
et al. I1996p . Hence, the a/Fe ratio quantifies the relative 
importance of Type II and Type la supernovae (Greggio 
& Renzini Matteucci & Greggio Pagel & 

Tautvaisiene ll995l Thomas et al. ll998p . and therefore car- 
ries information about the timescale over which star for- 
mation occurs. Thus, the a/Fe ratio can be considered to 
be an additional robust measure of the duration of the star 
formation process. Finally, a great deal of information on 
the formation process comes nowadays from simulations 
that are able to predict the radial distribution of the stellar 
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Fig. 19. Age, metallicity, and [a/Fe] gradients vs. the central velocity dispersion logcTc. The plotted stellar population 
parameters are derived from the combined analysis of the (H/3, <Fe>, Mgb), (R'-fp, <Fe>, Mgb), and (H^^, <Fe>, 
Mgb) index spaces (see text for details). Triangles are galaxies with relatively high density parameter (log(/9)>-0.4), 
and diamonds the other objects, while crosses mark the whole sample. The solid line is the linear best fit performed to 
all the galaxies. The dashed and dotted lines are the best fits to the galaxies located in relatively high and low density 
environments, respectively. 



populations within the galaxy (Larson ri974al Larson [T975l 
Carlberg [TOMl Bekki & Shioya ITMSl Chiosi & Carraro 
Kawata & Gibson HOOl Kobayashy HDUll). In the 
following we will summarize all the information we have 
deduced from the analysis of our sample, and see what sce- 
nario is more suited for early-type galaxies in low density 
environment. 

6.1. The scaling relations for field early-type galaxies 

An important result of our analysis is the derivation of cor- 
relations between the stellar population parameters and 
the galaxy potential well as a function of environment for 
our sample of galaxies. Other studies in the literature have 
presented scaling relations for early-type galaxies on the 
basis of different samples and/or stellar population mod- 
els, among which the most recent are those of Thomas et 
al. T05), Denicolo et al. (|2005b[ DOS), Gallazzi et 

al. (20011 G05) and Clemens et al. (l2006l C06). 

T05 have assembled three different data sets: G93, 
Beuing et al. ([2002]) and Mehlert et al. (f2000l [2003]) . Their 
sample contains objects located in both "high" and "low" 
surface density environments. It is worth noticing that 
H/? emission corrections have been applied by G93 and 
by Mehlert et al. ((20001 l2003l) but not by Beuing et al. 
(|2002p . T05 generated mock catalogues starting from the 



velocity dispersion distribution of the sample and assum- 
ing suitable scaling relations between the SSP parameters 
and the velocity dispersion (allowing for dispersion and for 
observational errors). The scaling relations are then mod- 
ified until the best "by-eye" fit to the indices (II/3, Mgb, 
and <Fe>) is obtained. 

DOS have derived ages, metallicities, and [a/Fe] ratios 
at re/8 for a sample of 83 early- type galaxies essentially 
in groups, the field, or isolated objects, which consist of 
S2 elliptical galaxies and 31 bulges of SOs or early- type 
spirals. Their results are based on the comparison of the 
measured II/3, Mgb, Mg2, <Fe> and [MgFe] indices with 
the TMB03 models. The 11(3 values used in their analysis 
were corrected for emission contamination through the Ha 
line (see also Sect. 2 of this paper). 

GOS presented stellar metallicities and ages for a mag- 
nitude limited sample of ~ 26,000 early-type galaxies 
drawn from the Sloan Digital Sky Survey Data Release 
Two (SDSS DR2). Their analysis rests on the recent pop- 
ulation synthesis models of Bruzual & Chariot p003p and 
is based on the simultaneous fit of five spectral absorption 
features (among which are II/3 and 116a + II7A), which de- 
pend only weakly on the a/Fe element abundance ratio. 
Before comparison with models, spectra were corrected for 
nebular emission lines. In their analysis [a/Fe] ratios have 
not been derived. 
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C06 have used the SSPs presented in this work to ana- 
lyze a volume limited subsample of about 4000 early-type 
galaxies in the Third Release of SDSS catalogue. By con- 
sidering 11 Lick indices at once, C06 have derived mean 
age, metallicity, and enhancement variations as a function 
of the dispersion velocity and environment. 

When comparing our results with the above literature 
we should keep in mind that: (a) the simple stellar popu- 
lation models adopted are different (in particular TMB03 
models are based on different stellar evolutionary tracks 
and on a different index dependence with element abun- 
dance with respect to our new models) ; (b) different line- 
strength indices are used in the analysis: in particular, 
while T05 and DG5 adopt a single age indicator, the H/3 
index, G05, and this work are based on the use of multiple 
age indicators (H/3, H(5, and H7), which allow us to mini- 
mize uncertainties due to emission contamination. Finally, 
COG use a simultaneous least-square fit to 11 index deriva- 
tives that takes their age/metallicity /enhancement (and 
also [C/Fe]) sensitivity into account; (c) Our sample is bi- 
ased not only towards early- type galaxies in LDE, but also 
towards galaxies whose spectra show emission lines, since 
we selected them to find ISM traces. T05 divide their sam- 
ple into low and high density classes, but in both works the 
classification is performed on the basis of a subjective es- 
timate of galaxy density (while we use density parameter, 
which is a properly measured volume density). The D05 
sample contains galaxies biased toward low density envi- 
ronment and G05 and C06 samples are drawn from the 
SDSS, which contains galaxies belonging to a wide range 
of environments, but predominantly composed of galax- 
ies in low density regions. Finally, DOS and G05 provide 
only qualitative trends of the derived stellar population 
parameters, while C06 provide differential variations as a 
function of velocity dispersion. In the following we will 
perform a detailed comparison of the results obtained in 
this work with those obtained in the above quoted pa- 
pers. However, for the sake of conciseness, we will limit 
the graphical comparison of the scaling relations to the 
work of T05, while we will only quote the results obtained 
by the other authors. 

6.1.1. The age scaling relation 

In the top panel of Fig. [20l we compare the age scaling 
relations obtained by T05 for the galaxies in low den- 
sity environments with our results for the total sample, 
for HDE galaxies and for LDE galaxies, respectively. It 
is worth noticing that the fit of T05 refers to a group 
of "passively" evolving galaxies (in the low density envi- 
ronment subsample) selected in the H/3-[MgFe] diagram, 
while we have selected, in our field sample, a subsample 
of "passively" evolving galaxies by considering a density 
threshold limit above which there is no evidence of reju- 
venation episodes (see Sect. 5.1.2). Thus, when comparing 
to T05, we consider their scaling relations for the galax- 
ies in a low density environment and our results for the 
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Fig. 20. Scaling relations obtained in this work compared 
with those of Thomas et al. (|2005p for their low density 
environment galaxies. Solid lines refer to Thomas et al. 
(12005p . while dashed, dot-dashed, and dotted lines refer 
to our HDE and LDE galaxies and to the total sample. 



HDE. The most noticeable difference between the age-CTc 
relation obtained in this work and that of T05 is in the 
slope. T05 find increasing ages with velocity dispersion 
(slope ~ 0.32), while we derive a flat relation (slope of ~ 
0.03). However, we notice that the linear correlation coef- 
ficient is poor in both cases (Eq. (f3T|) ). indicating that the 
spread of the data around the relation is large, and, given 
the limited number of objects, the error on the coefficients 
is also large. We also notice that if we consider the whole 
HDE + LDE sample, where young galaxies are present, 
we get a positive trend with a slope of 0.43. However, as 
already discussed in Sect. 5.1.1, a Spearman rank order 
test does not provide significant correlation between age 
and velocity dispersion. 

Our results are more consistent with those of G05 and 
C06. G05 do not find any significant variation of age with 
stellar mass (in their early-type sub-sample), except for 
a population of young metal-poor galaxies in the lowest 
mass bins. C06 find a significant flattening of the age above 
a ~ 170km/s. D05 estimate an average age of 5.8 Gyr for 
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Es and 3 Gyr for SOs, with hints of a trend for massive 
galaxies to be older. 

6.1.2. The metallicity scaling relation 

The relation between metallicity and velocity dispersion 
is shown in the central right panel of Fig. [201 The slope 
we have derived (~ 0.64) falls between the one derived 
by T05 (-0.57) and that derived by C06 (0.76). However, 
there is a significant zero point offset between T05 and our 
relation, with the metallicity derived by T05 being about 
0.15 dex larger than ours. At 100 km s^^, T05 find an 
average metallicity for galaxies in low density environment 
that is 26% larger than solar. Instead, we find that at the 
same velocity dispersion the average metallicity of old field 
galaxies should be —80% of the solar one. Our derived 
relation is in good agreement with metallicity estimates 
of stars of our own Bulge: the average metallicity of the 
Bulge is found to be [Fc/H]~-0.2 (Ramirez et al lMHH Mc 
William & Rich 1994j and, with its velocity dispersion of 
o-c =il04 km s^i (Blum. IT995j) . it faU just slightly below 
our relation. The zero point offset is most likely caused by 
the different evolutionary tracks adopted as fully discussed 
in Sect. 3. 

Apart from this offset, the estimate of the slope of the 
relation seems very robust, because it is quite indepen- 
dent from the adopted sample and from the analysis pro- 
cess. The linear correlation coefficient and the Spearman 
test indicate that the relation between metallicity and 
<jc is fairly tight. If we add the additional evidence that 
the slope seems to be almost independent from the envi- 
ronment (Eq. ([311) and T05), we may conclude that the 
metallicity-CTc relation is perhaps witnessing one of the 
tightest links between stellar populations and structural 
parameters. 

A metallicity increase with galaxy mass is found by 
D05 and G05 as well. However, the metallicities derived 
by D05 are larger when compared to our work and that of 
G05, likely because of the use of TMB03 models: at cr — 
300/cm/s DOS derive metallicities above 0.3 dex, while the 
G05 metallicities are peaked around logZ/Z© ~ 0.15 and 
our relation gives logZ/Z© ~ 0.22. We find only a slight 
dependence of the metallicity-cr relation on the richness 
parameter, in agreement with T05, while C06 find none, 
in their much larger sample. 

6.1.3. The a-enhancement scaling relation 

Our a-enhancement scaling relation is compared with the 
one obtained by T05 in the lower left panel of Fig. [20l 
Our derived slope (0.47) falls between that found by T05 
(0.28) and that of C06 (0.74) and agrees weU with that 
derived by D05 (0.59). Our [a/Fe] ratio coincides with 
that of T05 at crc=250 km s"\ while at CTc^IOO km s~i it 
is ~ 0.07 dex lower. The a-enhancement scaling relation 
is almost independent of the richness parameter, in very 
good agreement with C06 and T05. 



As discussed in Sect. 3.9, there are no significant dif- 
ferences in the [a/Fe] ratios predicted by TMB03 and our 
models. Moreover, since C06 have also used the models 
presented here, we may exclude that the different slopes 
found by T05, C06, and us are due to the models. The only 
noticeable differences are in the data samples and in the 
adopted methods. In particular, the use of different index 
sets coupled with the different sensitivity of the indices to 
the enhancement may explain the significant variation of 
the slope found by different authors. 

6.2. Gradients and galaxy formation 

Besides global stellar population parameters and their 
relation with galaxy structural parameters, population 
gradients within galaxies give one of the most strongest 
constraints on galaxy formation. Spectrophotometric and 
chemical properties at various locations within a galaxy 
are in fact closely related to the process of galaxy for- 
mation and evolution. Numerical simulations of dissipa- 
tive collapse of galaxies including star formation defi- 
nitely predict strong radial gradients in chemical enrich- 
ment (Larson IT974al[T975l Carlberg [19841). The predicted 
metallicity gradients are as steep as A log Z/ A log ~ - 
0.35 (Larson [T974al), ~ -1.0 (Larson [19751), and - -0.5 
(Carlberg ll984|) . A further prediction of monolithic mod- 
els is a steepening of the metallicity gradient with the po- 
tential well of the galaxy, which increases from about zero 
(in low mass galaxies) to ~ -0.5 in the most massive ones 
(Carlberg [HHH Bekki & Shioya[l999l Chiosi & Carraro 
[20021 Kawata & Gibson [20031 Kobayashi[2Q04l). 

On the contrary, the dissipationless collapse models 
predict no gradients in chemical enrichment (Gott 119731 
I1975P and the occurrence of major mergers predicts a sig- 
nificant dilution (Whiten[980! Bekki & ShioyalgM!) . More 
recently, Kobayashi (2004) has simulated the formation 
and chemodynamical evolution of elliptical galaxies in- 
cluding radiative cooling, star formation, feedback from 
Type II and la supernovae and stellar winds, and chem- 
ical enrichment. Galaxies are supposed to form through 
the successive merging of galaxies with various masses, 
which varies between a major merger at one extreme and 
a monolithic collapse of a slow-rotating gas cloud at the 
other extreme. They predict an average metallicity gra- 
dient of A log Z/ A log ~ — 0.3 ± 0.2 and no correlation 
between metallicity gradients and mass. The variety of 
gradients stems from the difference in the merging his- 
tories; in line with previous findings, galaxies that form 
monotonically have steeper gradients, while galaxies that 
undergo major mergers have shallower gradients. 

If metallicity gradients are more directly related to 
the degree at which star formation has reprocessed the 
gas inside the galaxy, [a/Fe] gradients give very strong 
constraints about the duration of the chemical enrich- 
ment process at different radii. According to Pipino & 
Matteucci (|2004|) . an outside- in formation scenario would 
predict [Mg/Fe] ratios increasing with galaxy radius (and 
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slightly older ages in the external regions compared to the 
central ones). On the other hand, an inside-out forma- 
tion, as suggested to explain the abundance gradients in 
the Milky Way (e.g., Matteucci & Francois 11^55)) . would 
produce a decrease of the [Mg/Fe] ratio with radius, since 
the outermost regions would evolve slower than the inner 
ones. 

From an observational aspect, evidence of metallicity 
gradients comes from the increase of line-strength indices 
(Carollo et al. 119931 Davies et al. [T9931 Sagha et al. [20001 
Trager et al. I2000al Mehlert et al. [MI51 Wu et al. [MIS]) 
and the reddening of the colors (e.g., Peletier I1990p to- 
wards the center of early-type galaxies. However, there is 
not yet a clear picture of the radial gradient behavior with 
respect to velocity dispersion. Some authors have shown 
that elliptical galaxies with larger values of the central 
Mg2 index tend to have steeper Mg2 gradients (Gorgas, 
Efstathiou & Aragon-Salamanca 119901 Carollo, Danziger 
& Buson FTOMl Gonzalez & Gorgas [lMg| . On the con- 
trary, Kobayashi & Arimoto (|1999p did not find any cor- 
relation between gradients and physical properties in the 
80 early-type galaxies of their study. Along the same line 
are the studies of Proctor & Sansom ^20021) and Mehlert 
et al. (pOOBj) . More recently, Forbes et al. ([2005) have re- 
analyzed the data of Coma cluster ellipticals (Mehlert et 
al. 120031 Sanchez-Blasauez l2004p and have found evidence 
for the presence of stronger metallicity gradients in more 
massive ellipticals. The only significant gradient shown by 
our study is that of the metallicity, while age and a en- 
hancement seem to be flat within Re/2, though with a sig- 
nificant dispersion fFigs. [TSl and [T9| . The average metal- 
licity gradient is AlogZ/Alog(r/re) ~ —0.21. There is 
also evidence for a trend of the Z gradient with central 
velocity dispersion, in the sense that more massive galax- 
ies tend to have shallower gradients (Fig. . We thus do 
not find evidence for a steepening of the metallicity gradi- 
ent toward larger masses, even if we consider only galaxies 
in the relatively richer environment. 

The flat distribution of a elements, on average, in spite 
of a well-established metallicity gradient, is a little intrigu- 
ing. The latter suggests that stellar populations are not 
mixed within re/2 and that metal enrichment was more 
efficient in the central region than in the outskirts of the 
galaxies. An outside- in formation scenario predicts the ex- 
istence of gradients in both metallicity and a enhance- 
ment, while merging would smear out stellar populations 
and dilute any gradient. Our observations indicate that 
the star formation proceeded on typical lifetimes not sig- 
nificantly different across re/2, but evidently with a larger 
efficiency in the center. It is easy to show, however, that 
the observed flat gradient in a enhancement is within the 
prediction of galaxy models. Using GALSYNTI^, we have 
run a chemical evolution model with parameters suited for 
a typical massive early-type galaxy (briefly, high SFR for 

2 GALSYNTH is a web based interface for 
the chemo-spectrophotometric code GRASIL at 
|http:/ /web. pd. astro. it/galsynth/index.php^ 



a fraction of a Gyr and passive evolution thereafter) and 
we have obtained that the run of [Mg/Fe] between [Z/H]=- 
0.5 and [Z/H]=0.5, where the majority of stars are formed, 
is well approximated by 

[Mg/Fe] = 0.46 - 0.13 x [Z/iJ], (34) 

Thus the expected gradient in [Mg/Fe] corresponding to 

AlogZ/Alog(r/re) 0.21 is A [Af g/Fe]/A log(r/re) ~ 

0.03 in very good agreement with what was observed. To 
find full support for an outside- in formation scenario, more 
external regions should be observed and analyzed. 

6.3. Rejuvenation 

The lack of environmental effects on the metallicity-iTc 
and [a/FeJ-CTc relations suggests that young ages are not 
the result of a more prolonged star formation in the low 
density environment. Moreover, NGC 3136, NGC3607, 
NGC7135, and NGC 6776 are massive galaxies, metal rich 
and a-enhanced, for which simple stellar population fits 
provide ages of only a few Gyrs. If recently formed and 
on the short timescale suggested by their a-enhancement, 
these objects would appear like the very powerful ultra lu- 
minous infrared galaxies (ULIRGs) observed in the local 
universe. However, it is easy to show that their luminos- 
ity would exceed by at least one order of magnitude that 
of the local powerful ULIRGs. In fact, recent determina- 
tions of star formation timcscales in these latter objects, 
by different methods, indicate that the gas content can- 
not sustain the observed luminosity for more than a few 
tens of Myr (e.g., Vega et al. [M351 Prouton et al. l^nUi|) . 
This implies SFR in excess of 10^ MQ /yr for the masses 
that correspond to the velocity dispersions of these galax- 
ies and, consequently, bolometric luminosities that by far 
exceed the typical ones in the local ULIRGs. 

All these considerations, together with the presence 
of morphology /kinematics signatures of recent interaction 
episodes in all the four quoted galaxies (see Appendix in 
Papers I and II) , strongly support the idea that the young 
ages obtained are due to recent rejuvenation episodes of 
star formation and do not correspond to their epoch of 
formation. If so, what is the mass involved in the rejuvena- 
tion episode? It is difficult to answer this question because 
the fading of the indices (H/3 for example) in a composite 
population depends on the fraction of the young compo- 
nent and the epoch of the event. We can only proceed 
with a statistical argument, though we are aware that the 
number of galaxies is not large. Looking to the data of 
Fig. [141 we determine that about 15% of all the objects 
look rejuvenated. The threshold limit for rejuvenation is 
taken from the age limit of the higher density objects. Let 
us suppose that, since the average epoch of formation, all 
the galaxies underwent a single rejuvenation event with 
a probability that follows the halo merging rate, namely 
--(l-fz)3-2 (Le Fevre et al. I2000p . The average epoch of 
formation is taken from the line in Fig. [14] indicating an 
average age of about 8 Gyr. 
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Assuming the above merging rate law, and a flat 
rim=0.3, r2A=0.7, and Hq — 70 km sec~^ Mpc~^ Universe, 
we determine that 15% of the total merging events since 
a look back time of 8 Gyr occur within the last 2.2 Gyr. 
This time must correspond to the time during which the 
H/3 index of a combined (old plus young) population re- 
mains above the selected threshold for rejuvenation and, 
of course, it will depend on the relative mass fractions of 
the old and young populations. Thus we have to seek for 
what fractional mass of the young component a combined 
population (8 Gyr plus the young component) maintains 
its H/3 index above the threshold value, for a duration of 
about 2.2 Gyr. 

After selecting a threshold value of H/3 ^^2.0 (from the 
data in Fig. [Ti|) . we find that the mass fraction of such a 
young population is ^ 25% of the total mass. We conclude 
that the data of Fig. [T3] are compatible with all galaxies 
having undergone a merger event since their formation, 
where not more than 25% of their total mass has been 
converted into stars. For a larger mass fraction, the frac- 
tion of rejuvenated objects we would see today would be 
larger. In other words, 75% of the total mass is definitely 
old, even if the galaxies suffered a major merger in their 
life. 

7. Summary and conclusions 

We have analyzed the stellar populations of the sample 
of 65 early-type galaxies mainly located in the field, pre- 
sented in Papers I and II. Our sample is biased towards 
early- type galaxies selected on the basis of ISM traces, in 
particular optical emission lines, tracing the warm ISM 
component. On the other hand, considering the incidence 
of galaxies showing a warm ISM in a randomly selected 
sample of early-type galaxies quoted by Falcon-Barroso et 
al. (2006), i.e., 75 %, (83% in the field), we are confident 
that our sample is not very dissimilar from those with 
which we compare our results. 

Stellar population parameters age, metallicity, and 
[a/Fe] ratio have been derived for all the apertures and 
gradients of the sample by the comparison of the Lick in- 
dices with new SSP models, presented in this work, which 
account for the presence of non-solar element abundance 
patterns. New SSP models have been developed for a wide 
range of ages (10^ — 16 x 10^ yr), metallicities (Z=0.0004, 
0.004, 0.008, 0.02, 0.05), and [a/Fe] ratios (0-0.8). The 
SSPs are based on the Padova stellar evolution tracks and 
accompanying isochrones (Bcrtelli et al. '1994; Bressan et 
al. [T994)) . which have been previously calibrated on color- 
magnitude diagrams of several star clusters in the Galaxy 
and in the Magellanic Clouds. 

The standard composition line-strength indices are 
computed on the basis of the fitting functions of Worthey 
et al. (lUM) and Worthey & Ottaviani (fTM7|) . The en- 
hanced SSPs are computed by applying an index correc- 
tion, which we derive following the main guidelines pro- 
vided by previous works in the literature (TMB03, Tantalo 
& Chiosi [^04ap . although with revised index dependence 



on the element abundance, to the standard SSPs. More 
specifically, we find that a logarithmic dependence of the 
index with element abundance provides a more satisfac- 
tory description of the index behavior. This evidence is 
also supported by preliminary results obtained running 
new model atmospheres and synthetic spectra with the 
code ATLAS12, that consistently computes both inte- 
rior and emergent spectra models with arbitrary abun- 
dances. The correction for the effects of a-enhancement is 
computed on the basis of the most recent specific index 
responses computed by Korn et al. (|2005p for different 
metallicities. 

The new a-enhanced SSPs have been implemented 
within an algorithm devised to derive stellar population 
parameters from the observed line-strength indices. The 
method, which rests on the probability density function, 
provides, together with the most probable solution in 
the (age, Z, [a/Fe]) space, the solutions along the age- 
metallicity degeneracy that are within 1 cr error from the 
observed index values. 

Summarizing our results: 

— The derived age distribution shows a wide spread, with 
SSP-equivalent ages from a few Gyr to 15 Gyr and with 
lenticular galaxies being generally younger than ellip- 
ticals (the average ages for the whole sample, E and SO 
are 8, 8.7, and 6.3 Gyr, respectively); the metallicity 
distribution shows a broad peak at < \Zl'H\ < 0.3; 
finally, the [a/Fe] ratio definitely presents a peak at 
^^0.22, with a narrower distribution than for age and 
metallicity. 

We also seek possible correlation of the stellar popula- 
tion parameters with both central velocity dispersion and 
local galaxy density. 

— Concerning the relation with the galaxy potential well, 
wc do not find the clear signature of a global trend of 
age with velocity dispersion; on the other hand, sig- 
nificant positive trends are derived for metallicity and 
a/Fe enhancement. 

— The robust metallicity-fTc relation testifies that chemi- 
cal enrichment is more efficient in more massive galax- 
ies and, within a universal initial mass function, the 
relation between the a/Fe enhancement and the ve- 
locity dispersion indicates that the overall duration of 
the star formation process is shorter within deeper po- 
tential wells. These two relations do not depend on 
galaxy morphological type, unlike the age-dc relation, 
indicating that the galaxy gravitational potential is the 
main driver of the chemical enrichment process of the 
galaxy. 

We analyzed the relations of the stellar population pa- 
rameters with the local galaxy density pxyz and we ob- 
tained the following results: 

— Concerning the age-pxyz relation, we derive a clear di- 
chotomy: while very young objects (from 1 Gyr to 4 



Annibali et al.: Nearby early-type galaxies with ionized gas. III. 



31 



Gyr, some of which have high CTc) are found in very low 
density environments {pxyz < 0.4), none of the galaxies 
in high density regions (40 % of the sample with mea- 
sured density) is younger than 4-5 Gyr. No trend of the 
chemical enrichment path is found with environmental 
conditions. The lack of environmental effect on the (a- 
enhancement)-tTc relation indicates that in very low 
environments rejuvenation episodes, rather than more 
prolonged star formation, are frequent. 

We analyzed the stellar population gradients within 
the single galaxies: 

— We definitely derive negative metallicity gradients 
(metallicity decreases from the central regions out- 
wards) of the order of A log 2'/A log(r/re) ^ —0.21, 
consistently with literature values. On average the a- 
enhancement remains quite flat within re/2, but this is 
consistent with the observed metallicity gradient. The 
only strong correlation is that of the metallicity gradi- 
ent with velocity dispersion: larger negative gradients 
are derived for less massive systems and become shal- 
lower at increasing galaxy mass. 

From the above results, we try to sketch a possible 
scenario for galaxy formation. The potential well of the 
galaxy seems to be the main driver of the chemical path 
during its formation, with more massive galaxies exhibit- 
ing the more efficient chemical enrichment and shorter star 
formation timescale. Furthermore, galaxies in denser envi- 
ronments are on average older than galaxies in less dense 
environments. This is in agreement with recent models of 
galaxy formation that combine feedback from QSO and 
SN, and cooling of the gas (Granato et al. [200T| [2004j) . 

Galaxies in less dense environments also show signa- 
tures of very recent rejuvenation episodes. By compar- 
ing the number of "young" objects with the total number 
of galaxies, and by means of simple two-SSP component 
models, we estimate that in these rejuvenation episodes 
(like major mergers), not more than 25% of the galaxy 
mass could be formed. Thus we estimate that, even if 
spheroidal galaxies are affected by major mergers during 
their evolution, about 75% of the galaxy mass is formed 
during the initial epoch of formation. 
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